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I.  INTRODUCTION 


Signal  processing  covers  a  wide  range  of  application  in 
such  diverse  fields  as  biomedical  engineering,  acoustics, 
sonar,  radar,  speech  communication,  data  communication,  and 
many  others.  As  an  example,  one  may  wish  to  remove  interfer¬ 
ence,  such  as  noise,  from  the  signal  or  to  modify  the  signal 
to  present  it  in  a  form  which  is  more  easily  interpreted  by 
an  expert.  As  another  example,  a  signal  transmitted  over  a 
communications  channel  is  generally  perturbed  in  a  variety  of 
ways,  including  channel  distortion,  fading,  and  the  insertion 
of  background  noise.  One  of  the  objectives  at  the  receiver 
is  to  compensate  for  these  perturbances.  In  each  case, 
processing  of  the  signal  is  required. 

Signal  proce-ssing  was  carried  out  typically  by  using 
analog  equipment  before  the  advent  of  digital  computers  which 
promoted  the  intensive  and  successful  development  of  digital 
signal  processing  in  the  past  decade.  At  present,  much  of 
the  digital  signal  processing  is  performed  on  general  purpose 
digital  computers  using  software.  However,  tremendous  progress 
has  also  been  made  in  other  types  of  signal  processors  using 
electronic  firmware  and  hardware.  Particularly,  using  elec¬ 
tronic  technologies,  a  new  class  of  sampled  analog  processors 
is  being  developed  using  charge  transport  devices.  Their 
signals  are  analog  but  the  independent  variables  are  discrete 
which  can  be  either  time  variable,  spatial  variables  or  other 
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transform  variables.  The  signal  processing  implementations 
are  illustrated  in  Table  I. I. 

Signal  processors  can  be  grouped  into  two  classes: 
filtering  and  spectral  analysis.  This  classification  is 
highlighted  in  Figure  1.1.  All  of  these  signal  processing 
systems  require  the  basic  essential  functions,  such  as  cor¬ 
relation,  convolution  and  transformation  (Fast  Fourier  trans¬ 
form).  In  other  words,  they  are  shown  to  be  centered  around 
four  basic  mathematical  operations:  delay  (or  shift),  mul¬ 
tiplication,  summation  and  generation  of  special  functions. 
These  operations  have  been  performed  by  digital  techniques 
and  it  was  probably  the  most  convenient  method.  However, 
the  new  sampled  analog  signal  processing  systems  only  require 
a  storage  for  discrete  time  analog  samples  and  not  discrete 
or  quantized  signal  amplitudes;  in  fact,  because  of  the 
computation  time,  frequency  selectivity,  and  part  cbunts,  it 
is  undesirable  to  quantize  the  signal.  Therefore,  discrete 
time  analog  devices  are  being  developed  to  perform  the  convo¬ 
lution,  correlation,  etc.  ,  without  the  need  to  convert  the 
signal  to  digital  format. 

There  are  several  distinct  advantages  of  using  the  discrete 
time  analog  devices  over  the  digital  techniques  to  perform  the 
same  algorithm: 

(1)  In  general,  the  discrete  time  device  system  can  perform 
the  same  computation  with  much  greater  speed  over  the 
digital  technique. 

(2)  The  device  stores  discrete  time  analog  amplitudes; 
therefore,  it  does  not  have  the  quantization  errors 
inherent  in  the  digital  system. 
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(3)  The  actual  multiplying  and  summing  operations  are 
easily  implemented,  and  with  fewer  components  than 
the  digital  system. 

To  illustrate  the  advantages  that  these  discrete  time 
devices  have  over  the  digital  technique,  a  simple  system 
shown  in  Figure  1-2  is  used.  This  is  a  model  of  a  simple 
single  pole,  first  order  discrete  time  filter.  To  implement 
this  system  with  a  discrete  time  device  such  as  the  SAD  (Serial 
Analog  Delay) ,  all  that  is  required  is  a  straight  substitution 
of  the  analog  device  for  the  delay  block,  an  operational 
amplifier  for  the  adder  and  a  potentiometer  proportional  to 
the  weighting  coefficient.  This  system  is  shown  in  Figure 
1-3. 

Comparing  this  system  implementation  to  the  digital  form 
in  Figure  1.4  demonstrates  the  difference  between  the  two 
systems.  It  requires  at  least  two  memory  registers,  one  to 
store  the  sampled  data,  the  other  to  store  the  multiplying 
constant.  In  the  figure  it  is  seen  that  more  devices  are 
required,  such  as  the  adder  and  the  costly  analog  to  digital 
converters,  as  well  as  the  digital  to  analog  converter  to 
complete  this  system. 

The  integrated  circuit  delay  devices  can  perform  in 
principle  most  of  the  necessary  algorithms  used  in  most  equip¬ 
ment  for  signal  processing  today.  However,  due  to  the  primi¬ 
tive  nature  of  these  new  classes  of  devices,  there  are 
limitations  in  some  fields  of  application. 

The  purpose  of  this  study  is  to  perform  the  evaluation  of 
Reticon  TAD-12  tapped  delay  lines  available  commercially  and 
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to  find  out  their  performances  and  limitations  in  order  to 
use  them  properly. 


TABLE  I . I  Signal  Processing  Implementation 


Signal 


Equipment 


Implementation 


Analog 


Digital 


Sampled 

Analog 


General  Purpose  Computer 
Special  Purpose  Computer 
Mini/Microcoraputer 
(Microprocessor) 

Signal  Processing  Function 

Modules : 

Multiplier 

Correlator 

Convolver  (matched 
filter) 

Filter 

Discrete  Fourier 
Transformer 

Frequency  synthesizer 

etc. 


Software 

Firmware 


Hardware 
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FIGURE  1.1  Classification  of  Signal  Processors 
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Time  Filter 


Input  >__vVvA. 


Multiplying 
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SAD  or 

1 

SAM 

1  Op.  Amp. 

<  Adder 

■>  Output 


FIGURE  1 . 3  Equivalent  Circuit  Implementation 


FIGURE  1 .4  Digital  Equivalent  Circuit  Implementation 
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II.  SAMPLED  ANALOG  NON- RECURSIVE  FILTER 


A.  PRINCIPLE 

The  transfer  function  of  a  digital  or  sampled  analog 
filter  can  be  expressed,  in  Z  domain,  by  one  of  the  follow¬ 
ing  equations. 


Recursive  filter:  H(Z)  =  - - 

1  j.  t'M  V  ■7~^ 

i  +  L-i  0^1 

1  m 

Non-Recursive  filter:  H(Z)  “  Zq  ^  ^ 


C2-1) 

C2-2) 


The  block  diagram  of  a  sampled  analog  recursive  filter 
is  illustrated  in  Figure  2-1:  _ 


E.  (z) 

in^  ' 


m 


FIGURE  2.1  Schematics  of  a  Recursive  Filter 
(Tapped  Delay  Line) 
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E 


in 


E  + 

out 

FIGURE  2,2  Schematic  of  a  Non-Recursive  Filter 
(Tapped  Delay  Line) 


In  the  case  of  causal  system  with  finite  duration  of 
impulse  response,  the  system  function  can  be  written  in  the 
form: 


H(Z) 


N--1 


_n 

h(n)  Z 


(2-3) 


{h(n)}  is  the  impulse  response  defined  over  the  time  interval 
0  £  nT  £  (N  -  1)T,  H(Z)  is  a  polynomial  in  Z  ^  of  degree  N  -  1. 
Thus,  H(Z)  has  N  -  1  zeroes  in  the  finite  Z-  plane. 

The  frequency  response  H(e^^)  is  obtained  by  substituting 
Z  by  e^^^. 
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(2-4) 


jwT.  N-1  -jwrtT 

H(e*'  )  =  Z  h(n)  e 

n=0 

We  recall  that  any  finite  duration  sequence  is  completely 
specified  by  N  samples  of  its  Fourier  Transform  [1] ,  so  that 
the  design  of  a  Finite  Impulse  Response  Filter  may  be 
accomplished  by  finding  either  its  impulse  response  coeffi¬ 
cients  or  N  samples  of  its  frequency  response.  Both  methods 
are  discussed  in  the  following  section. 


B.  LINEAR  PHASE  FINITE  IMPULSE  RESPONSE  FILTER 

In  many  applications,  e.g.,  speech  processing,  data 
transmission,  it  is  desirable  to  design  filters  to  have 
linear  phase.  In  this  way,  signals  in  the  passband  of  the 
filter  are  reproduced  exactly  at  the  filter  output  except 
for  a  delay  corresponding  to  the  slope  of  the  phase.  One 
of  the  most  important  features  of  FIR  systems  is  that  they 
can  be  designed  to  have  exactly  linear  phase.  The  unit 
impulse  response  for  a  causal  FIR  system  with  linear  phase 
has  the  property  that: 

h(n)  =  h(N  -  1  -  n)  (2-5) 


To  see  that  this  condition  implies  linear  phase,  we  write 
equation  2-4  as: 


H(2)  = 


CN/2)-l 

Z 

n=0 


(N/2)-l 

Z„ 


h(n)  Z'^  + 
h(n)  + 


N-1 

-z- 

n=N/  2 


(N/2)-l 

Z 

n=0 


h(n)  Z 

h(N-l-n) z" (N-1-^) 
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where  N  is  assumed  to  be  even.  Using  equation  2-5,  one  can 
write : 

.n  -fN-l-n'l 

H(Z)  =  I  h(n)[Z  n  +  2  (2-6) 

n=0 


If  N  is  odd,  it  is  shown  that 


Evaluating  equations  2-6  and  2-7  for  Z  =  e^'^,  one  obtains, 
for  N  even, 

H(e>)  =  ,-jw[(N-l)/2](W/’-l  2h(„)  cos  [w(n  - 

n=0  ^ 

(2-8) 

and  for  N  odd. 


[CN-3)/2]-l  , 

+  Z  2h(n)  cos  [W(n  -  ^V^)  ]  }  (2-9) 

n=0  ^ 


In  both  cases  the  sums  in  brackets  are  real,  implying  a 
linear  phase  shift  corresponding  to  delay  of  (N-l)/2  samples. 
For  N  odd,  the  phase  shift  corresponds  to  an  integer  number 
of  samples  delay.  For  even  N,  the  delay  is  an  integer  plus 
one-half  of  sampling  period.  This  distinction  between  odd 
and  even  values  of  N  is  often  of  considerable  importance  in 
design  and  realization  of  FIR  filter.  Examples  of  impulse 
response  having  linear  phase  are  shown  in  Figure  2-3. 
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FIGURE  2.3  Typical  Impulse  Response  for  Linear  Phase  FIR  Filters 

C.  DESIGN  TECHNIQUE  OF  FIR  FILTERS 

1 .  Determination  of  Impulse  Response 

Since  the  frequency  response  of  H(e^^)  of  any  digital 
or  sampled  analog  filter  is  periodic  in  frequency,  it  can  be 
expanded  in  Fourier  series  of  the  form: 

HCe^'^)  =  h(n)  e’^^”  (2-10) 

n=-“ 

where  the  sequence  h(n)  playing  the  role  of  the  ’’Fourier 
coefficients,”  i.e., 

h(n)  =  ^  HCe^"")  dw  (2-11) 
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One  can  easily  recognize  that  the  coefficients  of  the 
Fourier  series,  h(n) ,  are  identical  to  the  impulse  response 
of  a  digital  filter.  There  are  two  difficulties  with  the 
representation  of  (2-10)  for  designing  FIR  filters.  First, 
the  filter  impulse  response  is  infinite  in  duration  since  the 
summation  in  (2-10)  extends  to  ±0°  .  Second,  the  filter  is 
unrealizable  because  the  impulse  response  begins  at  -«>;  i.e., 
no  finite  amount  of  delay  can  make  the  impulse  response 
realizable . 

One  possible  way  of  obtaining  an  FIR  filter  that 
approximates  H(e^^)  would  be  to  truncate  the  infinite  Fourier 
series  (2-10),  at  n  =  ± M.  But  direct  truncation  of  the 

series  leads  to  the  well  known  Gibbs  phenomenon.  However, 
the  non-uniform  convergence  phenomenon  can  be  moderated 
through  the  use  of  a  less  abrupt  truncation  of  the  Fourier 
series  by  applying  the  so-called  "windowing”  or  "weighting 
function. " 

2 .  Effect  of  Window 

A  more  successful  way  of  obtaining  an  FIR  filter  is 
to  use  a  finite  weighting  sequence  W(n) ,  called  window,  to 
modify  the  Fourier  coefficients  h(n)  in  (2-10)  to  control  the 
convergence  of  the  Fourier  series.  The  technique  of  windowing 
is  illustrated  in  Figure  2-4.  At  the  top  of  this  figure  is 
shown  the  desired  periodic  frequency  response  H(e'^'^)  and  its 
Fourier  series  coefficients  {h(n)}.  The  next  row  shows  a 
finite  duration  weighting  sequence  W(n)  with  Fourier  trans¬ 
form  W(e^^)  .  1^(0^^^),  for  most  reasonable  windows,  consists 
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of  a  central  lobe  which  contains  most  of  the  energy  of  the 
windows  and  side  lobes  which  generally  decay  rapidly.  To 
produce  an  FIR  approximation  to  H(e^^),  the  sequence 


r.  _  j.h(n)  .  W(n) 
h(n)  =  {q 


-M  <  n  <  M 
elsewhere 


(2-12) 


is  formed.  The  third  row  of  Figure  2.4  shows  h(n)  and  its 
Fourier  transform  H(e^^) ,  which  is  the  convolution  of 
E:(e^^)  and  W(e^^)  ,  since  h(n)  is  the  product  of  the  sequences 
h(n)  and  w(n) .  The  last  row  of  Figure  2-4  shows  the  realiz¬ 
able  sequence  g(n) ,  which  is  a  shifted  version  of  h(n)  and 
may  be  used  as  the  desired  filter  impulse  response. 

As  seen  in  the  example  of  Figure  2-4,  there  are 
several  effects  on  the  resulting  frequency  response  of  win¬ 
dowing  the  Fourier  coefficients  of  the  filter.  A  major  effect 
is  that  discontinuities  in  HCe^'^)  become  a  transition  band 
between  values  on  either  side  of  the  discontinuity.  Since 
the  final  frequency  response  of  the  filter  is  the  circular 
convolution  of  the  ideal  frequency  response  with  the  window's 
frequency  response,  it  is  clear  that  the  width  of  these  tran¬ 
sition  bands  depend  on  the  width  of  the  main  lobe  of  W(e^'^)  . 

A  secondary  effect  of  windowing  is  that  ripple  from  the  side 
lobes  of  WCe^'^)  produces  approximation  errors  (ripple  in  the 
resulting  frequency  response)  for  all  w. 

There  are  many  windows  proposed  that  approximate  the 
desired  characteristics.  But,  in  general,  the  desirable 
window  should  have  the  following  characteristics: 
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FIGURE  2.4  Illustration  of  windowing 
Desired  periodic  frequency  response 
h(n)  Its  Fourier  series  coefficients 
w(n)  Finite  duration  weighting  sequence 
W(e^^)  Its  Fourier  transform 
h(n)  =  h(n)  .  w(n) 

H(e^’').  »  W(e^’) 
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1.  Small  width  of  main  lobe  of  the  frequency  response  of 
the  window  containing  as  much  of  the  total  energy  as 
possible . 

2.  Side  lobes  of  the  frequency  response  that  decrease  in 
energy  rapidly  as  w  tends  to  it. 

In  the  design  of  the  filters  in  this  study,  the 

Hamming  window  is  chosen.  It  is  specified  by  the  equation: 

,  0.54  -  0.46  cos  0  <  n  <  N-1 

wCn)  =  {  —  — 

0  ,  elsewhere  (2-13) 

Figure  2-5  illustrates  the  Hamming  window  and  the  rectangular 
window. 


W(n)  Rectan^lar 


FIGURE  2.5  Hamming  Window  and  Rectangular  Window 

3 .  Frequency  Sampling  Technique 

Figure  2-6  shows  an  arbitrary  frequency  response 
(solid  curve)  that  one  wants  to  approximate  and  a  sequence  of 
N  frequency  samples  H^(Figure  2-7)  that  can  be  represented  as 

H(k)  =  |H(k)|ej®^^\  k  =  0,  1...,  N-1 
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FIGIffiE  2.6  Desired  Continous  Frequency  Response 


FIGURE  2.7  Frequency  Samples 


Using  the  inverse  discrete  Fourier  transform,  a 
finite  duration  impulse  response  can  be  determined  from  HCk) 


h(n) 


H(k)  eH2’'/N)nk 

^  k=0 


(2-14) 


where  H(k)  can  also  be  written  in  the  form: 


H(k)  =  h(n) 
n=0 


(2-15) 


D.  APPLICATIONS 


1 .  Bandpass  Filter  Design  Algorithm 

In  the  analog  or  continuous  time  system,  the  bandpass 
filter  has  the  transfer  function  H(f). 


FIGURE  2.8  Bandpass  filter  frequency  response 

A,  f,  <  f  <  f. 

H(f)  =  {  ^  ^ 

0,  otherwise 

In  the  discrete  time  system,  a  rectangular  bandpass 
filter  has  the  frequency  response  H(e^^)  as  shown  in  Figure 
2-9.  It  is  periodic  of  period  2it  . 


0,  otherwise. 


Ws 

with  €2  -  sampling  frequency. 
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A 


PIGDEE  2.9  Discrete  Time  Bandpass  Filter  Frequency  Response 


Basing  upon  the  periodic  nature  of  a  discrete  time 
filter  frequency  response,  the  impulse  response  can  be  com¬ 
puted  by  finding  the  Fourier  coefficient,  using  equation  2-11. 

h(„)  .  f  A  eJ""  dw  (2-16) 

^TT  VT 

1 

1  ^2 

=  ±  f  A  (cos  wn  +  j  sin  wn)  dw. 
ir  w 

1 

Since  only  real  term  may  be  implemented,  the  series  coeffi¬ 
cients  are: 

Case  I :  n  ^  0 

1  ^2 

h(n)  =  i  A  cos  (wn)  dw 

=  —  (sin  Zirf  -  sin  Zirf  )  (2-17) 

irn  ^  2  1 

Case  2:  n  =  0 

Let  w  =  2it  (fo  + 

w  =  Ztt  (fo  + 

1  ^ 

One  can  write  using  equation  (2-16): 
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ur  ^  A  ,^^2  jwn  J 

h(n)  =  —  /  dw 

ir  Wj 

A  r  jw^n  jw,!!, 

"  “  -  eJ  .  ] 

”  T^T  (sin  ;rn4f) 

j  nTT  ^ 

=  2A  A  f (cos  2Trnfo  +  j  sin  2Tr  n  £o) 

(2-18) 

When  n  tends  to  0,  the  real  part  o£  equation  (2-18)  tends 
to: 

h(0)  =  2AA£(1.)C1.)  =  2AA  £  (2-19) 

Window  weighting  £unctions  are  used  to  inodi£y  the 
coe££icients  in  order  to  improve  the  £ilter  characteristics. 

I£  the  Hamming  window  is  used,  the  weighted  coe££icients 
A(n)  become 

A(n)  =  h(n)(0.54  +  0.46  cos  (^tt))  (2-20) 

j 

A  computer  program  shown  in  Appendix  C  was  established 
to  design  the  non-recursive  bandpass  £ilter.  It  was  based 
upon  the  equations  2-17,  2-19,  2-20. 

2 .  Prewhitening  and  Dewhitening  Filters  in 

Communication  System 

The  prewhitening  and  dewhitening  £ilters  are  also 
called  preemphasis  and  deemphasis  £ilters  in  the  communication 
£ield. 

Suppose  one  desires  to  transmit  a  baseband  signal 
using  FM  modulation  and  requires  the  best  possible  signal  to 
noise  ratio.  One  possible  solution  is  to  raise  the  level  o£ 
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modulating  baseband  signal  to  the  maximum  extent  possible 
in  order  to  modulate  the  carrier  as  vigorously  as  possible. 

But  the  modulating  signal  level  may  be  raised  only  until  the 
distortion  exceeds  a  specified  value  (maximum  allowable  value) . 

However,  the  baseband  signal  happens  to  be  an  audio 
signal;  it  turns  out  that  something  further  can  be  done.  An 
audio  signal  usually  has  the  characteristic  that  its  power 
spectral  density  is  relatively  high  in  the  low  frequency 
range  and  falls  off  rapidly  at  higher  frequencies.  For  ex¬ 
ample,  speech  has  little  power  spectral  density  above  about 
3  kHz.  Music  extends  farther  into  the  high  frequency  range, 
the  feature  still  persists  that  most  of  its  power  is  in  the 
low  frequency  region.  As  a  consequence,  when  one  examines 
the  spectrum  of  the  sidebands  associated  with  a  carrier  which 
is  frequency  modulated  by  an  audio  signal,  one  finds  that  the 
power  spectral  density  of  the  sidebands  is  greatest  near  the 
carrier  and  relatively  small  near  the  limits  of  the  allowable 
frequency  band  allocated  to  the  transmission.  The  manner  in 
which  one  may  take  advantage  of  these  spectral  features, 
which  are  characteristic  of  audio  signals,  in  order  to  improve 
the  performance  of  an  FM  system  is  shown  in  Figure  2-10. 

At  the  transmitting  end,  in  Figure  2-10,  the  baseband 
signal  m(t)  is  not  applied  directly  to  the  FM  modulator  but 
is  first  passed  through  a  filter  of  transfer  characteristic 
Hp (w) ,  so  that  the  modulating  signal  is  mp(t).  The  modifi¬ 
cation  introduced  into  the  baseband  signal  by  the  first 
filter  is  undone  by  the  receiver  filter  which  follows  the 
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discriminator  and  has  transfer  characteristic  1/Hp(w) .  The 
noise  passes  through  only  the  receiver  filter  which  may  be 
used  to  suppress  the  noise  to  some  extent. 

The  selection  of  the  transfer  characteristic  Hp(w) 
is  based  on  the  following  considerations.  At  the  output  of 
the  demodulator,  the  spectral  density  of  the  noise  increases 
with  the  square  of  the  frequency,  as  shown  in  Figure  2-11. 

It  is  given  by  the  equation: 

Gn(f)  “  ^  ■  |f  I  1  I  (2-21) 

where  a  =  constant  proportional  to  limited  amplitude 

of  the  carrier  obtained  with  the  hard  limiter 
in  the  demodulator 
A  =  amplitude  of  the  carrier 

n  =  power  spectral  density  of  the  white  noise 

B 

over  the  range  |f|  <  y  • 


FIGURE  2.11  Power  Spectral  Density  at  the  Output 
of  an  FM  Demodulator 
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Hence  the  receiver  filter  will  be  most  effective  in 
suppressing  noise  if  the  response  of  the  filter  falls  off 
with  increasing  frequency.  In  such  a  case  the  transmitter 
filter  must  exhibit  a  rising  response  with  increasing  fre¬ 
quency.  The  premodulation  filtering  in  the  transmitter,  to 
raise  the  power  spectral  density  of  the  baseband  signal  in 
its  upper  frequency  range,  is  called  preemphasis.  The  fil¬ 
tering  at  the  receiver  to  undo  the  signal  preemphasis  and 
to  suppress  noise  is  called  deemphasis. 

Referring  to  Figure  2-10,  one  requires  that  the  nor¬ 
malized  power  of  the  baseband  signal  m(t)  must  be  the  same 
as  the  normalized  power  of  the  preemphasized  signal  m^Ct). 

If  Gjjj(f)  is  the  power  spectral  density  of  m(t)  ,  the  density 
of  nip(t)  is  |Hp(f)l^  Gjjj(f)  ,  and  one  requires  that: 

fm  fm 

where 

fm  =  maximum  frequency  of  modulating  signal. 


In  the  absence  of  deemphasis,  output  noise  is: 


No  =  ^ 

With  deemphasis,  the  output  noise  is: 

fm 


f'od  ■  “’’'I  Ih-ttt 


df 


(2-23) 
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The  ratio,  R 


> 


R 


fmVS 


(2-24) 


Since  the  signal  itself  is  unaffected  in  the  overall 
process,  the  quantity  R  is  the  ratio  by  which  preemphasis- 
deemphasis  improves  the  signal  to  noise  ratio. 

The  prewhitening  filter  has  the  transfer  character¬ 
istic  illustrated  in  Figure  2-12.  The  dewhitening  filter 
has  the  transfer  characteristic  illustrated  in  Figure  2-13. 
The  relation  between  these  two  filters  is: 

HdCs^")  =  - (2-25) 

“  Hp(eJ’') 

where 

=  transfer  function  of  prewhitening  filter 
Hp(e^^)  =  transfer  function  of  dewhitening  filter. 

3 .  Prewhitening  Filter  Design  Algorithm 

a.  In  the  analog  or  continuous  time  system,  the 
prewhitening  filter  has  the  transfer  characteristic  shown 
in  Figure  2-12  and  has  the  transfer  function. 


H  (f) 
P 


R  ,  |£1  <  fo 

R  *  A(f  -  £o)  ,  fo  i  |£|  1  fp 

(2-26) 
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FIGURE  2.13 


f-  f  f  nonnalized 

0  p 

Dewhitening  Filter  Transfer  Function 


b.  In  the  sampled  analog  or  digital  system,  the 
frequency  response  of  the  prewhitening  filter  is  periodic. 
It  is  illustrated  in  Figure  2-14. 
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Analog  Prewhitening  Filter 


The  frequency  F  is  normalized  with  respect  to 
the  sampling  frequency  . 

The  design  of  a  linear  phase  prewhitening  filter 
consists  of  the  calculation  of  the  impulse  response  which 
can  be  computed  using  (2-11) 


h(n)  =  i  /  HCej")  dw 

i  II 

which  becomes,  using  ^  =  F 


h(n) 

h(n) 


2  R  dF  +  2/p^  [R  +  A(F-Fq)  ]  e^  ^’^”^dF 

^o 


U 


V 


(2-27) 


Since  h(n)  is  real,  one  can  compute  only  the  real 
part  of  U  and  V. 
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Case  1:  n  ^  0 


> 


F  =  ^ 


R 

Re[U]  =  2R  cos  (2TTnF)dF  =  —  sin  2TTn  Fq 


o 

F, 


Re[V]  =  2R  /  ^  cos  (2TTnF)dF  +  2A  /  ^  (F  -  Fq)  cos  (2TTnF)dF 


+  2A  (cos  2TTnFp  -  cos  2TTnFo) 

F_  Fr) 


AFq 

irn 


[sin  2TTnFo  -  sin  2iTnFp] 


h(n)  =  Re  [U]  +  Re  [V]  =  2^'2'j^2'  [cos  2TTn  Fp  -  cos  2iTn  Fq  ] 

(2-28) 

Case  2 :  n  =  0  ,  ~  ^s/2 


With  n  =  0,  (2-27)  becomes: 

Fp  Fp 

h(o)  =  2/  RdF  +  2/  [R  +  A(F  -  Fq)  ]dF 
op 
^0 

hCo)  =  [R  *  ^4:  (Fp  -  Fo)=l 


2Fp 


2Fp  '-"P 


(2-29) 


Taking  into  account  the  constraint  of  the  linear 
phase  FIR  filter: 

h(n)  =  h  (N  -  1  -  n) 

and  using  (2-28)  and  (2-29),  a  computer  program 
is  prepared  to  calculate  the  impulse  response  of 
the  filter,  a  Hamming  window  function  is  used  to 
improve  the  amplitude  and  phase  response  of  the 
prewhitening  filter. 
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4 .  Dewhitening  Filter  Design  Algorithm 

In  the  design  of  the  complementary  dewhitening  filter, 
one  can  use  the  frequency  sampling  technique.  Using  (2-15) 
and  the  given  impulse  response  of  the  prewhitening  filter, 
the  frequency  samples 

{Hp(k) }  ,  k  =  0,  1 _ ,  N-1 

can  be  computed  and  stored  in  the  memory  of  the  computer. 

Then  the  frequency  samples  of  dewhitening  filter 
Hjj(k)  can  be  computed  as: 

^D ~  H  Ck)'  *  k=0,  1,....,N-1 

Using  the  inverse  discrete  Fourier  transform  given 
by  (2-14)  ,  the  finite  impulse  response  of  the  dewhitening 
filter  can  be  obtained. 

Using  an  appropriate  window,  such  as  Hamming  window, 
the  frequency  response  of  the  dewhitening  filter  can  be 
improved  also. 
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III.  EXPERIMENTAL  EVALUATION  OF  TAD- 12 


A.  SYSTEM  CONFIGURATION 

Reticon  Corporation  markets  a  Tapped  Analog  Delay  Device 
which  provides  12  output  taps  with  the  following  time  delay 
configuration: 


*  PIGDHE  5*1  '  BqvilTalent  Circuit  of  TAD- 12 

The  first  delay  section  is  delayed  by  one  clock  period, 
the  rest  have  delay  time  of  two  clock  periods. 

The  delay  lines  are  structurally  organized  as  shown  in 
Figure  3-2.  One  delay  line  contains  N  capacitive  storage 
elements  with  two  I/O  access  switches.  One  access  switch 
connects  the  storage  element  to  a  common  input  line  and  the 
second  switch  connects  it  to  the  common  output  line.  One 
register  and  its  associated  multiplexing  switch  sequentially 
time  samples  the  analog  input  signal  onto  the  ”N"  capacitive 
storage  elements.  Information  is  read  onto  the  nth  element 

1“  H 

while  simultaneously  reading  out  the  n^ 


+  1  element. 


CLOCK  INPUT 


FIGURE  5»2  Organization  of  a  Single  Delay  Line 
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Each  element  has  an  output  buffer  amplifier  to  provide 
a  zero  order  hold  of  the  sampled  signal.  This  produces  a 
sampled  and  held  output. 

The  test  circuit  for  TAD-12  is  shown  in  Figure  3-3.  It 
consists  of  a  pulse  shaper  circuitry  and  a  DMOS  flip-flop. 
Because  the  tap  outputs  are  at  +5  volts  bias  level,  due  to 
the  TAD-12  circuitry,  the  resistor  2  kilo-ohms  of  the  oper¬ 
ational  amplifier  non- inverting  input  is  connected  to  a  +5 
volts  rather  than  ground.  The  output  of  the  summer,  there¬ 
fore,  rides  on  a  5  volts  bias  level. 

B.  FREQUENCY  DOMAIN  EVALUATION 

Using  the  circuit  board  shown  in  Figure  3-3,  the  frequency 
dependance,  the  non-uniformity,  the  harmonic  distortion  are 
evaluated  at  different  sampling  frequencies  for  three  differ¬ 
ent  tap  resistances: 

=  10  Kfi  ,  100  Kfi  ,  1  Mfi  . 

1 .  Frequency  Roll-off 

The  frequency  characteristic  of  the  sampled  and  hold 
process  is  of  the  ^  shape  with  x  =  ^  . 

X  15 

At  sampling  frequency  f  =  500  KHz,  for  R  =  10  Kfi, 

^  L. 

the  roll-off  is  close  to  the  theory.  When  the  tapping 
resistances  R^.^p  is  increased,  the  roll-off  is  more  serious 
and  wriggles  around. 

Figure  3-4  illustrates  the  magnitude  of  the  frequency 
response  of  Tap  #1,  Tap  #6  and  Tap  #11  outputs  using  10  kilo- 
ohms  resistor.  The  solid  curve  represents  the  theoretical 


response . 
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FIGURE  3.3  Test  Circuit  for  Reticon  TAD-12  Evaluation 
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+5V  -15V  +15V 


From  this  experiment  it  is  found  that  the  TAD- 12 

output  has  rather  complicated  frequency  variations  although 

s  in  X 

they  follow  the  theoretical  -  relation  roughly  if  the 

tapping  resistor  is  10  Kilo- ohms  and  the  sampling  frequency 
less  than  2.5  MHz.  The  loading  effect  is  much  more  pronounced 
if  the  tapping  resistors  are  equal  or  greater  than  100  Kilo- 
ohms  . 

2 .  Non-Uniformity  (Spatial  or  Fixed  Pattern  Noise) 

The  non-uniformity  of  different  taps  outputs  varies 
depending  on  the  tap  resistance,  Rj^,  the  sampling  frequency, 
f s ,  and  the  input  signal  frequency,  f.  The  general  result 
at  sampling  frequency  of  500  KHz  is  the  following: 

-  For  f  <  10  KHz,  the  non-uniformity  is  within  10% 

-  For  f  >  10  KHz,  the  non-uniformity  can  deteriorate 
to  15%  -  50%  range. 

Figures  3-10  through  3-17  illustrate  the  non-uniformity 
result  of  sampling  frequency  of  500  KHz,  for  different  input 
signal  frequencies  and  for  different  tap  resistors: 

=  1  Kj^,  10  Kn,  100  Kn,  1  Mn. 

At  high  sampling  frequencies,  using  tap  resistances 
less  than  100  Kn,  the  uniformity  of  tap  output  can  be  improved 
at  higher  input  signal  frequencies. 

At  fg  =  0.5  MHz,  Rj^  =  10  KJ^  and  input  signal  frequencies: 

-  f  <  50  KHz,  the  non-uniformity  is  within  10% 

-  f  >  50  KHz,  the  non-uniformity  varies  from  13%  to  25%. 
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At  £5  =  2.5  MHz,  Rj^  =  100  Kfi  and  input  signal  frequencies: 

-  £  <  80  KHz,  the  non- uniformity  is  within  101 

-  £  >  80  KHz,  the  non-uniformity  varies  from  11%  to  70% 

At  £g  =  2.5  MHz,  =  1  Mfi,  and  input  signal  frequencies: 

-  f  <  30  KHz,  the  non-uniformity  is  within  20% 

-  £  >  30  KHz,  the  non-uniformity  is  greater  than  30%. 

Figures  3-18  through  3-28  illustrate  the  non-uniformity 
data  at  sampling  frequency  equal  to  2.5  MHz,  for  different 
typical  tap  resistors  and  for  different  input  signal  fre¬ 
quencies  . 

Figures  3-24  through  3-40  illustrate  the  pictorial  pres¬ 
entation  of  Reticon  TAD- 12  performance  at  different  sampling 
f  re*quencies ,  using  different  tap  resistances. 

3 .  Temporal  Noise  Behavior  of  Different  Taps 

The  temporal  noise  of  different  taps  is  studied  at 
the  sampling  frequencies  higher  or  equal  to  2.5  MHz.  It 
varies,  depending  on  the  sampling  frequencies  and  the  tap 
resistances  used.  However,  it  should  be  pointed  out  that 
the  tap  #1  output  has  little  noise  at  all  sampling  frequen¬ 
cies  and  tapping  resistors  used.  It  is  shown  in  Figures  3-24, 
3-25,  3-26,  3-28,  3-40.  The  general  temporal  noise  behavior 
of  other  taps  is  shown  in  Table  III. I  and  Figures  3-24,  3-25, 

3-27,  3-35,  3-38. 


50 


TABLE  III. I  Noisiest  Tap  Outputs 


K 

tap 

1  Meg-ohm  100  Kilo-ohms  10  Kilo-ohms 

2.5  MHz  Tap  #10  Tap  #11  Tap  #12 

3.5  MHz  Tap  #10  Tap  #4  Tap  #2 

The  temporal  noise  of  different  taps  is  not  serious 
at  low  sampling  rate  (less  or  equal  to  1.25  MHz).  Figure 
3-41  illustrates  the  output  of  tap  #1  and  #12  at  sampling 
frequency  at  fg  =  1.25  MHz  using  1  Meg-ohm  tapping  resistors. 

4 .  Harmonic  Distortion 

The  second  harmonic  distortion,  data  is  measured  with 
the  input  frequency  10  KHz  for  different  input  amplitudes  at 
different  sampling  frequencies: 

fs  =  500  KHz,  1.25  MHz,  2.5  MHz 
and  using  different  tap  resistors: 

Rtap  =  10  Kfi,  100  1  MJ^. 

The  curves  for  the  quantity  20  log  y  fundam^^al" 
are  plotted  and  shown  in  Figures  3-42,  3-43,  3-44.  The  general 
results  of  the  above  measurement  are  summarized  in  Table  3. II. 

The  loading  effect  on  the  dynamic  range  of  the  TAD-12 
is  very  pronounced  at  high  sampling  frequency  (2.5  MHz).  In 
order  to  get  the  second  harmonic  amplitude  less  than  1%  of 
the  fundamental,  one  has  to  use  less  than  .12  volt  rms  input 
voltage  if  using  1  Mn  tap  resistor  and  less  than  2  volts  rms 
input  voltage  if  using  10  Kn  tap  resistor. 
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FIGURE  3.4  Amplitude  Frequency  Response  of  Reticon  TAD-12  Tapped  Delay  Line. 
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FIGURE  3*5  Amplitude  Frequency  Response  of  Reticon  TAD  -  12  Tapped  Delay  Line 
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FIGURE  5»6  Frequency  Response  of  Reticon  TAD-12  Tapped  Delay  Line  Showing 
Its  Bandwith,  Nonuniformities  Among  Taps  and  Loading  Effect 
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FIGURE  3»7  Amplitude  Frequency  Response  of  Reticon  TAI)-12  Tapped  Delay  Line  « 
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FIGURE  3.8  Amplitude  Frequency  Response  of  Reticon  TAD-12  Tapped  Delay  Line. 
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FIGURE  5.9  Amplitude  Frequency  Response  of  Reticon  TAD-12  Tapped  Delay  Line, 
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FIGURE  3*11  Non  Uniformity  Data,  f  =*  500  KHz 
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PIGDEE  5.12  Non  Uniformity  Data,  f^„„  -  500  KHz 
"tap"  ^inpat’ 


Tap  Output  Amplitude 

PIGDEE  3.13  Non  Uniformity  Data.  f„^„  -  500  KHz 

ScUiip 

E.  =10  Kilo-ohms,  100  KHz 

tap  input 
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FIGURE  3.14  Non  Ifciformity  Data.  f„^  “  500  KHz 
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FIGURE  3.15  Non  Uniformity  Data. 
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FIGURE  5*17  Non  Uniformity  Data,  f  =  500  KHz 
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Tap  Output  Amplitude 

FIGURE  5 •18  Non  Uniformity  Data,  f  *  2.5  MHz 

samp 

R.  -  10  Kilo-ohms.  f.  ,  50  KHz 

tap  input 
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FIGURE  3.19  Non  Uniformity  Data,  f  _  =  2.5  MHz. 
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Tap  Output  Amplitude 

FIGURE  5 *20  Non  Uniformity  Data,  f  =2.5  MHz 

^  samp 
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tap  input 


FIGURE  5.21  Non  Uniformity  Data, 

®tap“  Kilo-ohms,  ^input"* 
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^ap“  ■'  MeS-ohm. 


R.  =1  Meg-ohm.  f.  .=■  100  KHz 
tap  input 
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Input  Frequency  30  KHz  100  KHz 
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FIGURE  5.24  Pictorial  Presentation  of  TAD-12  Performance*  ^gajnp“  2. 5  MHz.  ^  Meg-ohm 

V  ,*  .06  volt  peak  to  peak.  Scale;  H  ®  10  microsec. /div;  V  «  .05  volt/div 
input 


20  KHz  50  KHz  100  KHz  I50  KHz 


A  A 


^tap”  Kilo-ohms;  ^^nput**  Volt  peak  to  peak(sinusoidal  input)  . 
Scale:  H*10  raicroaec./div, ,  V«,2  V/div, 


Tap  #1  Output 


Tap  #5  Output 


Tap  #6  Output 


Tap  #11  Output 


FIGURE  3*26  Pictorial  Presentation  of  TAI>-12  Performance. 

f  “3  MHz;  R.  =100  Kilo-ohms, 
samp  tap 

Vinput*  volt  peak  to  peak;  ^^nput^* 

Scale:  H  =  10  microsec. /div,  V  =  .2  volt/div 
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Tap  #1  Output 


Tap  #4  Output 


Tap  #2  Output 


Tap  #3  Output 


FIGURE  3.27 


Pictorial  Presentation  of  Reticon  TAD-12 


Performance,  f  =  3*5  MHz. 
samp 

7  .6  Volt  peak  to  peak, 

input 

Scale:  H  =  10  microsec. /div; 


-100  Kilo-ohms 

tap 


^ input 
V  =•  .2 


*  50  KHz. 
Volt/div 
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Tap  #1  Output 


Tap  #2  Output 


Tap  #3  Output 


Tap  #12  Output 


FIGURE  5»28  Pictorial  Presentation  of  Reticon  TAJ)-12, 


Performance 


V,  .6 

input 

Scale :  H  = 


.  f  -  4  MHz.  R,  - 
samp  tap 

V  peak  to  peak; 
r  10  microaec./div;  V  = 


100  Kohma. 
50  KHz. 

.2  V/div. 
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Tap  #1  Output 


Tap  #9  Output 


FIGURE  5.29 


Reticon  TAI)-12  Performance. 

f  =  4.5  MHz.  R.  =  100  Kilo-ohms, 
samp  tap 

Vinput”*  ^  peak  to  peak.  ^input°' 

Scale:  H  =  10  micro sec./div;  V  *  .2  ?/div. 


Tap  #1  Output 


Tap  ffO.  Output 


FIGURE  3 *50  Reticon  TAD- 12  Performance. 

f  =5  MHz.  R,  =100  Kilo-ohms, 
samp  tap 

^input“  ^  ^input’  50  KHz. 

Scale;  H  =  10  microsec./div;  V  =  .2  V/div. 
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Input 


Tap  #1  Output 


FIGURE  3«31  Reticon  TAI)-12  Performance. 

f  “2.5  MHz.  R.  “  10  Kilo-ohms 
samp  tap 

V-  .6  V  peak  to  peak.  f.  .=* 

input  ^  input 

Upper  Scale:  H  =  5  micro  sec. /div; 
Lower  Scale:  H  =  5  microsec. /div; 


FIGURE  3.32  Reticon  TAD-12  Performeince 

f  =2.5  MHz.  R,  =10  Kilo-ohms, 
samp  tap 

V.  .=  .6  V  pesUc  to  T)eak.  f.  .= 
input  ^  ^  input 

Upper  Scale:  H  =  .5  microsec. /div; 
Lower  Scale:  H  =  .5  microsec. /div; 


100  KHz. 

V  =  .2  V/div 

V  =  .1  V/div 


MHz. 

=  .2  V/div 
=  .1  V/div 
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Input 


Tap  #6  Output 


FIGURE  5»33  Reticon  TAI)-12  Performance. 

f  *  2.5  MHz.  R.  =10  Kilo-ohms, 
samp  tap 

Vinput*  ^  peak  to  peak,  ^^nput"  KHz. 
Upper  Scale:  H  =  1  microsec. /div;  V  =  .2  V/div 
Lower  Scale:  H  =  1  microsec. /div;  V  =  .1  V/div 


Input 


Tap  #6  Output 


FIGURE  3 ‘34  Reticon  TAI>-12  Performance. 

f  =2.5  MHz;  R.  =10  Kilo-ohms, 
samp  tap 

V.  .»  .6  V  peak  to  peaik.  f,  .=  1  MHz 
input  ^  ^  input 

Upper  Scale:  H  =  .5  microsec. /div;  V  =  .2  ?/div 
Lower  Seale:  H  =  .5  microsec. /div;  V  =  .1  v/div 
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Input 


Tap  #12  Output 


FIGURE  3.35 


Reticon  TAD-12  Performance. 


samp 


2.5  MHz.  R, 


tap 


10  Kilo-ohms. 


V.  .=  .6  V  peak  to  peaUc.  f. 
input  ^  ^  input 

Upper  Scale:  H  =  5  microsec. /div;  V 
Lower  Scale:  H  =•  5  microsec. /div;  V 


100  KHz 
»  .2  V/div 
-  .1  V/div 


Input 


Tap  #12  Output 


FIGURE  3.36  Reticon  TAD- 12  Performance. 

^samp'  ^-5  ■‘tap’  Kilo-ohms. 

’'input"  ''  300  KHs 

Upper  Scale:  H  *  2  microsec. /div;  V  »  .2  V/div 
Lower  Scale:  H  =  2  microsec. /div;  v  =  .1  V/div 
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Input 


Tap  #3  Output 


FIGURE  3*37  Reticon  TAD- 12  Performance. 


f  »  3  MHz.  R,  =10  Kilo-ohms, 
samp  tap 

V.  ,=>  .6  V  peak  to  peak;  f.  , 
input  ^  ^  input 

Upper  Scale:  H  =  2  microsec./div; 
Lower  Scale:  H  =  2  microsec./div; 


=  200  KHz 

V  =  .5  V/div 

V  =  .2  V/div 


Input 


Tap  #2  Output 


FIGURE  3*30  Reticon  TAD- 12  Performance. 

f  *  3* 5  MHz;  R,  =10  Kilo-ohms, 
samp  tap 

V.  ,=  .6  V  peak  to  peak;  f.  ,=  200  KHz 
input  ^  ^  input 

Scale:  H  =  2  microsec./div;  7  =  .5  V/div 
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Input 


Tap  #2  Output 


FIGURE  3*39  Reticon  TAD- 12  Performance. 

f  »  5  MHz.  R.  =10  Kilo-ohms 
samp  tap 

V,  .=  .6  V  peak  to  peak.  f.  .=  200  KHz 
input  ^  ^  input 

Upper  Scale;  H  =  5  micro sec./div;  7  =  .2  V/div 
Lower  Scale;  H  »  5  microsec. /div;  V  =  .01  V/div 


Input 


Tap  #1  Output 


FIGURE  3*40  Reticon  TAD-12  Performance. 

f  =5  MHz.  R.  =10  Kilo-ohms 
samp  tap 

7.  .=  .6  7  peak  to  peak;  f.  .=  200  KHz 

input  input 

Upper  Scale;  H  =  2  microsec. /div;  7  «  .2  7/div 
Lower  Scale;  H  =  2  microsec. /div;  7  =  .1  7/div 
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Input  Frequency  50  KHz  100  KHz 
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77 


20  log(V  2nd  harm*/  V  fund.) 
■  10  Kilo-ohms 
-100  Kilo-ohms 
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500  KHz  £c  =  1.25  NHz  =  2.5  MHz 
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c.  TIME  DOMAIN  EVALUATION 

The  measurement  procedure  is  to  use  a  periodic  input 
pulse.  Its  period  is  at  least  twice  the  total  delay.  The 
width  of  the  pulse  must  satisfy  the  condition: 

Y  of  clock  period  £  width  of  pulse  <  1  clock  period. 

The  pulse  amplitude  should  be  within  the  dynamic  range.  At 
the  output  of  the  summer  operational  amplifier,  one  can  ob¬ 
serve  twelve  output  amplitudes.  An  analog  sample  and  hold 
circuit  is  used  to  measure  each  output  amplitude. 

Figures  3-45  through  3-47  illustrate  the  impulse  response 
of  TAD- 12  performance  at  different  sampling  frequencies: 

fs  =  50  KHz,  1.25  MHz,  2.5  MHz 

and  tap  resistances: 

=  10  Kfi,  100  KQ,  1  Mn 

connected  to  tap  #1,  5,  6. 

1 .  Non- Uniformity 

The  non-uniformity  data  of  tap  outputs  for  the  impulse 
input  is  shown  in  Figures  3-48,  3-49,  using  10  Ki^  tap  resist¬ 
ances  at  different  sampling  frequencies: 

fs  =  25  KHz,  50  KHz,  500  KHz,  1.25  MHz. 

The  variation  of  the  tap  output  amplitudes  is  the 
following : 


4.21, 

at 

fs  -  25 

KHz 

1.8%, 

at 

£3  =  50 

KHz 

4.1%, 

at 

£3  =  500 

KHZ 

10.47% 

,  at 

£5  =  1. 

25  MHz. 
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Figures  3-50,  3-51  illustrate  the  output  of  12  taps 
using  =  10  sampling  frequencies: 

fg  =  50  KHz,  1.25  MHz. 

At  lower  sampling  frequency,  up  to  f^  =  1.25  MHz,  the 
non-uniformity  of  tap  output  is  less  than  10%;  when  the  samp¬ 
ling  is  increased,  the  uniformity  is  deteriorated  seriously. 
At  f^  =  2.5  MHz,  the  non-uniformity  of  different  taps  is 
higher  than  20%. 

The  interaction  between  adjacent  taps  can  also  affect 
the  tap  uniformity  at  high  sampling  frequencies.  Figure 
3-45  (c) ,  (d) ,  (ej ,  (f)  shows  these  effects.  For  example, 
in  Figure  3-45  (dj  at  f^  =  1.25  MHz,  the  tap  #1  output  can 
have  an  effect  on  tap  #2  before  its  amplitude  reaches  the 
zero  value.  Another  example,  in  Figure  3-45  (f) ,  at  f^  of 
2.5  MHz,  the  effect  of  tap  #1  output  can  extend  to  tap  #2 
and  tap  #3  before  its  amplitude  goes  below  the  fixed  pattern 
noise . 

2 .  Temporal  Noise 

The  noise  of  tap  outputs  can  be  grouped  into  two 
categories : 

-  fixed  pattern  noise 

-  random  noise. 

The  fixed  pattern  noise  is  periodic  with  a  period 
related  to  the  sampling  frequency  used.  At  f^  of  50  KHz, 
this  period  is  about  6.67  microseconds.  At  f^  of  500  KHz, 
it  is  about  80  microseconds,  but  the  number  of  periods  is 
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constant  between  tap  outputs.  For  example,  there  are  three 
periods  of  fixed  pattern  noise  between  tap  #1  output  and 
tap  #5  output.  This  is  illustrated  in  Figures  3-45  (a), 

(c) ,  3-46  (a),  (c)  ,  3-47  (a),  (c)  .  The  peak  to  peak  ampli¬ 
tude  of  this  noise  is  about  40  mV  for  the  10  Kfl  load  resistors 
and  about  25  mV  for  the  100  Kfl  load. 

The  random  noise  varies  from  one  tap  to  the  other. 

It  depends  also  on  the  clock  frequency  and  tap  resistors. 

The  random  noise  affects  seriously  the  tap  output,  except 
tap  #1,  at  high  tap  resistor  value  and  high  clock  frequency. 
For  example,  using  10  Kfi  resistors,  the  amplitude  of  random 
noise  is  40  mV  at  clock  frequency  of  f^  =  2.5  MHz  and  is 
100  mV  at  f^  =  5  MHz.  This  is  shown  in  Figure  3-45  (d) ,  (f) . 
The  tap  #5  output  is  deteriorated  badly  at  f^  =  5  MHz  for 
R  =  100  Kfi  and  R  =  1  Mf2  in  Figures  3-46  (d)  and 
3-47  (d) .  However,  the  noise  on  tap  #1  is  very  small  at  any 
condition.  This  is  illustrated  in  Figures  3-45  (d) ,  (f ) , 

3-46  Cd),  (f). 

The  jitter  amplitude  observed  at  the  top  of  tap  #1 
output  in  Figure  3-52  equals  to  80  microvolts  for  a  signal 
frequency  of  2.83  KHz,  f  =  100  KHz  and  for  R  =  10  Kfi. 

C  L 
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Figure  3-5  shows  the  magnitude  frequency  response 
of  Tap  #1,  Tap  #6  and  Tap  #8  outputs  using  100  Kilo-ohms 
resistors  and  1  Meg-ohm  resistors.  The  roll-off  starts  at 
about  20  KHz  input  frequency  for  the  =  1  Meg-ohm  curve 

and  at  about  60  KHz  input  frequency  for  the  R^^p  =  100  Kilo- 
ohms  curves . 

When  the  sampling  frequency  is  increased,  the  fre¬ 
quency  dependence  curves  are  not  improved  in  comparison  with 
the  theoretical  curve.  Figures  3-6  and  3-7  illustrate  the 
magnitude  frequency  response  of  different  tap  outputs  at 
sampling  frequency  equal  to  2.5  MHz  using  tap  resistors: 

R^^p  =  10  Kilo-ohms,  100  Kilo-ohms,  1  Meg-ohm. 

With  R^^p  =  10  Kilo-ohms,  the  frequency  roll-off  is 
about  150  KHz.  When  the  signal  frequency,  f,  equals  to  200 
KHz,  the  measured  curve  is  1  db  below  the  theoretical  curve. 
At  f  =  300  KHz,  it  is  3  db  below  the  theoretical  curve.  For 
the  R.^  =  1  Meg-ohm  curves,  the  frequency  roll-off  begins 

about  30  KHz. 

With  R^^p  =  100  Kilo-ohms,  the  response  curve  begins 
to  roll  off  at  the  input  frequency  f  =  80  KHz.  At  f  =  100 
KHz,  the  response  curve  is  2  db  below  the  theoretical  curve. 

Figures  3-8  and  3-9  illustrate  respectively  the 
magnitude  frequency  response  of  TAD- 12  output  using  10  Kilo- 
ohms  tap  resistors  at  sampling  frequency  f^  equal  to  2.75  MHz 
and  using  1  Meg-ohm  tap  resistor  at  fg  equal  to  1.25  MHz. 
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50  KHz  1.25  MHz  2.5  MHz 
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Lower  Scale:  H=1  microsec./div.  V=>.02  V/div.;  (d)  Scale:  H-1  micros ec./div.  V-.1  V/div. 

(e)  Upper  Scale:  H=2  microsec./div.  V=1  V/div. ,  Lower  Scale:  H=2  microsec./div.  V-.02  V/div. 

(f)  Upper  Scale:  H=10  microsec./div.  V-.05  V/div,,  Lower  Scale:  H=10  microsec./div.  V».1V/div 
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(b)  Scale:  H»10  mlcrosec./div,  V-.5  V/div. ;  (c)  Upper  Scale:  H-1  microsec. /div.  V-1  V/div. , 

Lower  Scale:  H=>1  microsec, /div.  V»,05  V/div.;  (d)  Scale:  H»5  microsec, /div,  V=,02  V/div,  | 

(e)  Upper  Scale:  H=2  microsec. /div,  V”2  V/div, ,  Lower  Scale:  H«2  microsec, /div,  V*.01  V/div. ; 

(f)  Scale:  H»2  microsec. /div.  V»,1  V/div, 
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(b)  Scale:  H*50  micro sec •/div.  V».1  V/div. ;  (c)  Upper  Scale:  H“.5  microsec. /div.  V»2  V/div. , 
Lower  Scale:  H-.5  microsec./div.  V-.002  V/div. ;  (d)  Scale;  H»5  microsec. /div.  V«.01  v/div.; 

(e)  Upper  Scale;  H=.2  microsec. /div.  V=1  V/div. ,  Lower  Scale;  H=*.2  microsec./div.  V’=.002  V/div. » 

(f)  Upper  Scale;  H-10  microsec./div.  V-.01  v/div. ,  Lower  Scale;  H-10  microsec./div.  V=.02  V/div. 


PIGT3HE  3»48  Non  Dhifonnity  Data.  Input  Pulse:  Amplitude»1  Volt, 
Width=.2  microsec..  Tapping  Resistor=10  Kilo-ohms 

(a)  P  =  25  KHz 

'  ^  samp 

(b)  P  =  50  KHz 

'  samp 


88 


Tap  output  amplitude 


FIGURE  3 *49  Non  Uniformity  Data*  Tapping  Resistora^lO  Kilo-ohms 

(a)  Input  Pulse r  Amplitude=1 . 1  Volt,  Width=,6  microsec. 


F^^^  -500  KHz 
samp 

(b)  Input  Pulse;  Amplitude-1.1  Volt, 


Width-. 5  raicrosec. 


F 

samp 


1.25  Mz 


8.9 


Output 


Input 


?IGDRE 


Output 


Input 


FIGURE  5 


.50  Impulse  Response  of  Reticon  TAD- 12  showing 
12  Tap  Outputs,  ^-tap*  Kilo-ohms. 

Sampling  Frequency  =  50  KHz. 

Upper  Scale:  H  =  50  microsec. /div;  V«  .5  V/div 
Lower  Scale;  H  =  50  microsec. /div;  7=1  V/div 


,51  Impulse  Response  of  Reticon  TAD-12  Showing 
12  Tap  Outputs,  Kilo-ohms. 

Sampling  Frequency  =  1.25  MHz. 

Upper  Scale;  H  =  2  microsec. /div;  V  =  .2  V/div 
Lower  Scale;  H  =  2  micro sec. div/;  7=1  7/div 
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Output 


Tap  #1  Jittering 
Waveform 


FIGURE  5»52  Upper  Trace:  Impulse  Response  of  12  Tap  Outputs 

Scale:  H  »  .1  msec./div;  V  =  ,05  V/div 
Lower  Trace:  Jittering  Waveform  of  Tap  #1 

Scale:  H  =  ,1  msec./div;  V  »  ,02  Y/div 
Sampling  Frequency  »  50  KHz 
Tapping  Resisteuice  *10  Kilo-ohms, 
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IV.  INVESTIGATION  OF  FILTER  I  -  PREWHITENING  FILTER 


A.  THEORY 

Based  on  the  design  algorithm  described  in  II.D.l,  a 
computer  program  shown  in  Appendix  A  is  written  to  give  the 
impulse  response  of  the  prewhitening  filter  for  the  required 
frequency  characteristics.  Using  this  program,  a  prewhitening 
filter  is  designed  for  the  following  conditions: 

f^  =  2.1  KHz  f  =  1/2  f  =  625  KHz 

o  p  s 

R  =  0.1  Slope  A  =  1 

Length  of  filter:  M  =  11  taps 


The  tapping  coefficients  improved  by  Hamming  window  are 
listed  in  Table  IV. I. 

The  magnitude  frequency  response  of  the  theoretical  pre¬ 
whitening  filter  is  shown  in  Figure  4-1.  It  is  plotted  in 
linear  scale  on  Y-axis  versus  normalized  frequencies, 
f/fg  on  X-axis. 
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TABLE  IV. I  Theoretical  Tapping  Coefficients  §  Resistors 


k 

-5 

-4 

-3 

-2 

-1 

0 

1 

2 

3 

4 

5 


-0.00091 

0.0 

-0.01219 

0.0 

-0.21407 

1. 

-0.21427 

0.0 

-0.01219 

0.0 

-0.00091 


(Kilo-ohms) 

-16478.516 

open 

-1225.517 

open 

-65.005 

10. 

-65.005 

open 

-1225.517 

open 

-16478.516 


A,  Measured 
k 

0.02564 
0.15385 
0.05128 
0.  0 

-0.23076 

1. 

-0.19231 
0.02564 
0.02564 
0. 02566 
-0.07692 


B.  EXPERIMENT 


1 .  Tapping  Resistors  According  to  Theoretical 

Design  with  ^min  -  10  Kilo-ohms 

A  minimum  value  of  10  Kilo-ohms  is  set  for  the  tapping 


resistances.  Appendix  F  indicates  the  following  expression 
for  Rj^, 


R,  =  (lOK  +  r)  0 


^max 


)  -  r 


where  r  is  the  internal  output  resistance.  For  the  TAD-12, 
r  is  about  5  Kilo-ohms.  The  tapping  resistors  corresponding 
to  the  designed  prewhitening  filter  are  given  in  Table  IV. 1 
and  used  to  implement  the  filter  using  the  11  taps  of 
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Reticon  TAD-12.  Since  their  theoretical  coefficients 
for  k  higher  than  5  have  negligible  values,  the  implementa¬ 
tion  of  the  filter  on  a  longer  delay  line,  number  of  taps 
higher  than  11,  is  not  necessary. 

The  measured  impulse  response  of  this  filter  is 
illustrated  in  Figure  4-2  at  three  different  sampling  fre¬ 
quencies  : 

fg  =  50  KHz,  500  KHz,  1.25  MHz 

The  pulse  inputs  used  have  the  amplitude  and  width  listed 
in  Table  IV. II. 

TABLE  IV. II  Pulse  Inputs 
fg  Amplitude  Width 

50  KHz  1.5V  8  micros ec. 

500  KHz  1.5  V  .8  microsec. 

1.25  MHz  1.5  V  .3  microsec. 

Using  the  sample  and  hold  circuit,  the  impulse 
response,  A^^  is  measured  and  the  result  is  listed  in  Table 
IV.  1.  The  Aj^  measured  and  the  theoretical  Aj^  are  quite 
different  for  some  taps.  The  maximum  deviation  occurs  be¬ 
tween  tap  #7  and  #8  outputs.  Its  peak  amplitude  is  about 
301  of  the  output  at  tap  #6. 

The  theoretical  magnitude  frequency  response  and 
measured  response  curves  for  three  different  sampling  fre¬ 
quencies  are  plotted  on  Figure  4-3.  The  response  curves  for 
f^  of  1.25  MHz  and  500  KHz  show  some  wiggling  shape  due 
probably  to  the  non-uniformity  of  tap  output  and  the  clock 
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noise  which  cause  the  unsymmetrical  tapping  coefficients. 

The  response  curve  for  f  equal  to  50  KHz  is  rather  satis¬ 
factory  if  it  is  multiplied  by  a  factor  of  1.07. 

2 .  Tapping  Resistors  Adjusted  to  Yield 

Better  Impulse  Response.  Rmin  =  10  Kilo-ohms 

Since  only  three  taps  have  strong  contribution  to 

the  filter,  the  adjustment  is  done  only  on  the  taps  #5,  #6, 

#7  by  using  the  sample  and  hold  circuit.  The  other  taps  are 

left  open.  The  tapping  resistors  used  are  listed  in  Table 

IV. III. 


TABLE  IV. I II.  Tapping  Resistors 

k  A.  (Kilo-ohms) 

^  fs  =  50  KHz  fs  =  500  KHz  fs  =  1.25  MHz 


1 

-0.21427 

28 

50 

52 

0 

1. 

10 

10 

10 

1 

-0.21427 

28 

35 

100 

The  amplitude  and  width  of  inputs  pulse  are  listed  in 
Table  IV. II. 

The  impulse  response  of  this  filter  is  illustrated  in 
Figure  4-4,  at  three  different  sampling  frequencies. 

The  theoretical  and  measured  curves  of  prewhitening 
filter  magnitude  frequency  response  are  shown  in  Figure  4-5. 
The  response  curves  for  f^  equal  to  50  KHz,  500  KHz  and  1.25  MHz 
still  show  some  wiggling.  The  gain  of  these  curves  are  lower 
than  the  last  curves  in  Figure  4-3.  If  these  three  curves 
are  multiplied  by  a  factor  of  1.22,  they  are  somewhat  satis¬ 
factory  in  comparison  with  the  theoretical  curve. 
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3 .  Tapping  Resistors  Adjusted  to  Yield 

Better  Impulse  Response.  Rmin  =  50  ohms 

A  minimum  value  of  50  ohms  is  set  for  the  tap  #6. 

The  taps  #5  and  #7  are  adjusted  by  using  sample  and  hold 

circuit,  to  yield  better  impulse  response.  The  other  taps 

are  left  open.  The  tapping  resistors  used  are  listed  in 

Table  IV. IV. 

TABLE  IV. IV.  Tapping  Resistors 

k  A]j  (Kilo-ohms) 

fg  =  50  KHz  fg  =  500  KHz  fg  =  1.25  MHz 


1 

-0.21427 

00 

2.8 

00 

0 

1. 

0.05 

0.05 

0.05 

1 

-0.21427 

2.8 

1.8 

11.5 

The  pulse  inputs  used  have  the  amplitudes  and  widths 
listed  in  Table  IV. II. 

Figure  4-6  illustrates  the  impulse  response  of  the 
filter  at  sampling  frequencies; 

fg  =  50  KHz,  500  KHz,  1.25  MHz 

The  theoretical  and  measured  curves  of  prewhitening 
filter  magnitude  frequency  response  are  shown  in  Figure  4-7. 
The  wiggling  of  the  three  measured  curves  is  much  reduced  in 
comparison  with  the  previous  curves.  If  the  experimental 
curve  f2=l*25  MHz  is  multiplied  by  a  factor  of  1.2  and  curves 
fg  =  50  KHz,  500  KHz  multiplied  by  a  factor  of  1.35,  the 
agreement  is  somewhat  satisfactory. 
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4 .  Tapping  Resistors  Obtained  by  Experimental 

Cut  §  Try  Based  on  Frequency  Response  Adjustment 

The  procedure  consists  of  the  use  of  a  chirped  wave¬ 
form  from  DC  to  1  MHz  as  the  input  of  the  filter,  and  adjust 
the  displayed  frequency  response  to  get  the  minimum  wiggling 
by  adjusting  the  tapping  resistances.  The  result  obtained 
is  shown  in  Figure  4-9.  The  wiggling  of  the  three  response 
curves  is  much  reduced.  Raising  these  curves  by  a  factor  of 
1.27  can  improve  the  agreement  between  the  theoretical  and 
measured  curves. 

The  impulse  response  corresponding  to  these  curves 
is  unsymmetrical  and  has  the  values  listed  in  Table  IV. V. 

TABLE  IV. V  Tapping  Coefficients 
k  Aj.  Rk  (Kilo-ohms) 

-1  -0.1694  24.812 

0  1.  .05 

1  -1  .05 

The  amplitude  and  width  of  pulse  inputs  used  are 
listed  in  Table  IV. II. 


97 


COL  dCZ  OO'i  006  OUi  CIO  010  OW 


FIGURE  4.1  Theoretical  Amplitude  Frequency  Response  of  Prewhitening 
Filter.  Length  =*  11.  Hamming  window. 

X  -  Scale  “  .1  units/inch  for  normalized  frequency 

Y. -  Scale  »  .2  units/inch  for  Amplitude  frequency  response. 
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(c)  Upper  Scale:  H-.3  miorosec./div.  V»1,5  V/div.;  Lower  Scale:  H»5  microsec. /div.  V»2  v/di 


H(f)4 
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FIGURE  4.3  Prewhitening  Filter  Magnitude  Frequency  Response. 

Tapping  Resistors  According  to  Theoretical  Design,  Kilo-ohms 


50  KHz  500  KHz  1.25  MHz 
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(c)  Upper  Scale:  H»5  microsec./div.  V'*,02  V/div.;  Lower  Scale:  H*5  mlcroseo./div,  V«1  V/div 
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FIGURE  4/. 5  Prewhitening  Filter  Magnitude  Frequency  Response. 

Tapping  Resistors  Adjusted  to  Yield  Better  Impulse  Response.  R^^^  -  10  Kilo-ohms 
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(c)  Upper  Scale:  H=*5  microaec*/div*  V«.5  V/div.;  Lower  Scale:  H»5  laioroaec./div.  V®2  V/div 
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FIGURE  ^-7  Prewhitening  Filter  Magnitude  Frequency  Response. 
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(c)  Scalet  H*»5  mlcrosec./div.  V*»1  V/div, 


H<f) 


106 


FIGURE  i$.9  Prewhitening  Filter  Magnitude  Frequency  Response. 

Tapping  Resistors  Obtained  by  Cut  and  Try  Based  on  Frequency  Domain  Response. 


V.  INVESTIGATION  OF  FILTER  II  -  DEWHITENING  FILTER 


A.  THEORY 

The  design  algorithm  of  the  dewhitening  filter  is 
described  in  II. D. 2.  The  computer  program  in  Appendix  B  is 
used  to  calculate  the  impulse  response  of  the  filter  from 
the  given  frequency  response.  A  dewhitening  filter  is 
designed,  based  on  the  impulse  response  of  the  prewhitening 
filter  of  paragraph  IV. A.  The  frequency  characteristics 
are  the  following: 

fo  =  2.1  KHz  "  I  =  625  KHz 

Length  of  filter  M  =  11  taps 


frt  f  f  normalized 

0  p 

The  tapping  coefficients  are  improved  by  a  Hamming  window 
and  by  a  factor  of  1.13  and  are  listed  in  Table  V.I. 

The  magnitude  frequency  response  of  the  theoretical  pre¬ 
whitening  filter  is  shown  in  Figure  5.1.  It  is  plotted  in 
linear  scale  on  Y-axis  versus  normalized  frequencies, 
f/fs  on  X-axis. 
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TABLE  V . I .  Theoretical  Tapping  Coefficients  §  Resistors 


k  theor.  (Kilo-ohms)  (Kilo-ohms) 

with  Rjnin  “  K  ohms,  with  Rmin  =  50  ohms 


-5 

0.00113 

13269.336 

4464.026 

-4 

0.00632 

2368.417 

794.050 

-3 

0.01156 

1292.577 

431.851 

-2 

0.04846 

304.661 

99.252 

-1 

0.19681 

71.215 

20,659 

0 

1. 

10. 

.  05 

1 

0.19681 

71.215 

20.659 

2 

0.04844 

304.661 

99.252 

3 

0.01156 

1292.577 

431.851 

4 

0.00632 

2368.417 

794.050 

5 

0.00113 

13269.336 

4464.026 

EXPERIMENT 

1. 

Tapping  Resistors 

According  to 

Theoretical 

Design  with  Rmin 

=  10  Kilo-ohms 

A  minimum  value  of  10  Kilo-ohms  is  set  for  the 
tapping  resistors.  Using  the  expression  of  Rj^  given  in 
Appendix  F,  the  tapping  resistances  are  calculated  and 
listed  in  Table  V.I  in  the  third  column.  They  are  used  to 
implement  the  dewhitening  filter  using  11  taps  of  Reticon 
TAD- 12. 

The  impulse  response  of  this  filter  is  illustrated 
in  Figure  5-2  at  three  different  sampling  frequencies: 

fs  =  50  KHz,  500  KHz,  1.25  MHz. 
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TABLE  V.II  Theoretical  ^  Measured  Coefficients 

\ 

k  A]^  theor.  A]^  Measured 

fs  =  50  KHz  fs  =  500  KHz  fg  =  1. 25  MHz 


-5 

0.00113 

0.16216 

0.10869 

0.11667 

-4 

0.00632 

0.16216 

0.17391 

0.20830 

-3 

0.01156 

0.16216 

0.15217 

0.25 

-2 

0.04844 

0.2027 

0.19565 

0.2708 

-1 

0.19681 

0.39119 

0.41304 

0.41667 

0 

1. 

1. 

1. 

1. 

1 

0.19681 

0.39189 

0.39130 

0.41667 

2 

0.04844 

0.18918 

0.13043 

0.0 

3 

0.01156 

0.16216 

0.13043 

0.16667 

4 

0.00632 

0.17567 

0.13043 

0.16667 

5 

0.00113 

0.16216 

0.10869 

0.16667 

The  impulse 

inputs  are  noted  in  Table 

IV. II  for 

different 

sampling 

frequencies . 

The  tapping 

coefficients 

are  measured  with  the  sample  and 

hold  circuit 

and  the  result 

are  listed 

.  in  Table 

V.II.  The  theoretical  and  measured  Aj^ 

are  quite 

different 

for  some  taps 

due  to  the 

fixed  pattern 

noise  and 

non-uniformity  of  tap  outputs.  The 

maximum  noise 

amplitude 

is  observed  before  the 

tap  #1  output  and  has 

amplitude  of  about  30%  of  Aq. 

The  theoretical  magnitude  frequency  response  and 
measured  response  for  three  different  sampling  frequencies 
are  plotted  in  Figure  5-3.  The  response  curves  for  f^ 
of  500  KHz  and  1.25  MHz  are  somewhat  in  agreement  with  the 
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theoretical  curve.  Their  difference  is  about  12%  for  f 
less  than  0.2  f^  and  it  becomes  worse  than  40%  for  input 
signal  frequency  higher  than  0.25  f^.  For  f^  of  50  KHz, 
the  result  is  quite  satisfactory. 

2 .  Tapping  Resistors  Adjusted  to  Yield  Better 

Impulse  Response.  Rmin  =  10  Kilo-ohms 

Since  the  coefficients,  for  k  greater  than  3, 
have  negligible  values  in  comparison  with  the  noise,  the 
adjustment  is  done  only  on  seven  tap  outputs.  The  adjust¬ 
ment  is  quite  delicate  for  the  tap  outputs  which  are  greater 
than  or  equal  to  5,  because  of  the  interaction  between  them. 
Experimental  tap  resistances  are  listed  in  Table  V.III.  A 
slight  modification  has  to  be  done  on  the  values  of  some 
tap  resistors  in  order  to  get  better  impulse  responses  for 
different  sampling  frequencies.  The  sample  and  hold  circuit 
is  used  to  make  the  adjustment. 

The  amplitude  and  width  of  inputs  pulse  used  are 
listed  in  Table  IV. II. 

The  impulse  response  of  this  filter  is  illustrated 
in  Figure  5-4  at  three  different  sampling  frequencies.  The 
maximum  noise  occurs  near  tap  #1  output;  its  peak  amplitude 
is  about  19%  of  A^. 

The  theoretical  and  measured  curves  of  dewhitening 

filter  magnitude  frequency  response  are  shown  in  Figure  5-5. 

The  wiggling  of  the  response  curves  is  reduced  in  comparison 

with  the  curves  in  Figure  5-3,  but  the  difference  between  the 

theoretical  and  measured  curves  (f  =1.25  MHz  and  500  KHz) 

'■  s 
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is  still  high,  greater  than  30%,  for  the  input  signal 
frequency  f  higher  than  0.25  f ^ .  At  low  sampling  frequency, 
50  KHz,  the  result  is  quite  satisfactory. 

TABLE  V.III.  Adjusted  Coefficients  §  Resistors 


k  Ak  Rk 

f  =50  KHz  f  =500  KHz  f  =1.25MHz  f  =50  KHz  f^  =  500  KHz  f  =1 . 25MHz 
s  s  s  s  s  s 


-5 

0.09836 

0.02083 

0.125 

13269. 

13269. 

13269. 

-4 

0.10655 

0.02083 

0.1 

2368. 

2368. 

2368. 

-3 

0.10655 

0.04167 

0.075 

1292.5 

1292.5 

1292.5 

-2 

0.09836 

0.0625 

0.075 

300. 

300. 

300. 

-1 

0.18852 

0.20833 

0.175 

110. 

109. 

110. 

0 

1. 

1. 

1. 

5. 

9.5 

5. 

1 

0.18852 

0.20833 

0.2 

110 

110. 

110. 

2 

0.09836 

0.0625 

0.1 

311. 

311. 

320. 

3 

0.10655 

0.02083 

0.075 

1292.5 

1295.5 

1292.5 

4 

0.11675 

0.02083 

0.075 

2368. 

2368. 

2368. 

5 

0.10655 

0.02083 

0.075 

13269. 

13269. 

13269. 

3.  Tapping 

Resistors  According 

to  Theoretical  Design. 

Rmin 

50  ohms. 

A  minimum  value 

of  50  ohms 

is  set  for 

the  tap  #6. 

Using 

the 

formula  described  in 

Appendix  F 

,  one  can 

calculate  the 

theoretical  tapping  resistances  whose  values  are  listed  in  the 
fourth  column  of  Table  V.I. 

The  input  pulses  have  the  amplitude  and  width  listed 
in  Table  IV. II.  Figure  5-6  illustrates  the  impulse  response 


111 


of  the  filter  at  sampling  frequencies: 


fg  =  50  KHz,  500  KHz,  1.25  MHz 
The  measured  coefficients  are  listed  in  Table  V.IV: 


TABLE 

V.IV.  Measured 

Coefficients  (R„,- 
^  min 

=50  ohms) 

k 

f  =50  KHz 
s 

^  -^k 

fg=500  KHz 

f^=1.25  MHz 
s 

-5 

0.11111 

0.06667 

0.22826 

-4 

0.11111 

0.06667 

0.08695 

-3 

0.11111 

0.06667 

0.09782 

-2 

0.16667 

0.13333 

0.15217 

-1 

0.27778 

0.26667 

0.34782 

0 

1. 

1. 

1. 

1 

0.22223 

0.26667 

0.48913 

2 

0.13888 

0.13333 

0.20652 

3 

0.11111 

0.06667 

0.13043 

4 

0.11111 

0.06667 

0.04347 

5 

0.11111 

0.06667 

0.07608 

The  maximum  peak  amplitude  of  noise  measured  is  about 
131  of  for  f^  equal  to  1.25  MHz  and  500  KHz  and  19%  of  A^ 
for  f^  equal  to  100  KHz. 

The  theoretical  and  measured  curves  of  dewhitening 
filter  magnitude  frequency  response  are  shown  in  Figure  5-7. 
The  response  curve  for  f^  equal  to  50  KHz  is  similar  to  the 
theoretical  one.  The  curves  f^  equal  to  1.25  MHz  and  500  KHz 
have  some  wiggling  effect  for  input  frequencies  higher  than 


.29  f3. 
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000  000  OiO  OiS  020  02b  030 


FIGURE  5.1  Amplitude  Frequency  Response  of  Dewhitening  Filter. 
Length  *  11  /  Hamming  window 

X  -  Scale  =  .1  units/inch  for  normalized  frequency 
Y  -  Scale  *  .5  units/inch  for  amplitude  frequency 
response. 
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(c)  Upper  Scale;  H»5  micro sec. /div.  V«,01  V/div. ;  Lower  Scale:  H-5  microsec. /div,  V«1  V/di 
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FIGURE  5.3  Dewhitening  Filter  Frequency  Response. 

Tapping  Resistors  According  to  Design.  “10  Kilo-ohms 


50  KHz  5D0  KHz  1,25  MHz 
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(c)  Upper  Scale:  H-5  microaec./div,  V-.01  V/div, ;  Lower  Scale:  H»5  micro sec. /div.  V»2  V/di 
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FIGURE  5.5  Dewhitening  Filter  Magnitude  Frequency  Response. 

Tapping  Resistors  Adjusted  to  Yield  Better  Impulse  Response  .  -  10  Kilo-ohms 


50  KHz  500  KHz  1.25  MHz 
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(c)  Upper  Scale:  H=*5  microsec. /dlv,  V».02  V/div.;  Lower  Scale:  H-5  microsec. /div.  V-1  V/di 
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VI.  INVESTIGATION  OF  FILTER  III  -  BANDPASS  FILTER 


A.  THEORY 

Based  upon  the  design  algorithm  described  in  II. D. 3,  a 
computer  program  shown  in  Appendix  C  is  written  to  calculate 
the  impulse  response  of  the  bandpass  filter  based  on  the 
required  frequency  response.  Using  this  program,  a  bandpass 
filter  is  designed  for  the  following  conditions: 

=  .1996  £2  =  .2004 

Center  frequency  f ^  =  . 2  Q  =  250 

Amplitude  A  =  1. 

Number  of  taps  M  =  11. 


The  tapping  coefficients  improved  by  Hamming  window  are 
listed  in  Table  VI. I. 

The  magnitude  frequency  response  of  the  theoretical  band¬ 
pass  filter  is  shown  in  Figure  6-1,  in  linear  scale,  and  in 
Figure  6-2,  in  logarithmic  scale  versus  normalized  frequencies, 

f/fs- 
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TABLE  VI. I.  Theoretical  Coefficients  and  Resistors 


k 

=  k)  Rv(R 
k  min  k^  min 

-5 

0.09863 

147.083 

46.201 

-4 

0.07378 

198.307 

63.447 

-3 

-0.38389 

34.074 

8.155 

-2 

-0.59148 

20.360 

3.538 

-1 

0.28647 

47.362 

12.628 

0 

1. 

10. 

.05 

1 

0.28647 

47.362 

12.628 

2 

-0.59148 

20.360 

3.  538 

3 

-0.38389 

34.074 

8.155 

4 

0.07378 

198.307 

63.447 

5 

0.09863 

147.083 

46.201 

EXPERIMENT 

1. 

Tapping  Resistors 

According  to  Theoretical 

Design  with  Rmin  ' 

=  10  K  ohms. 

A  minimum  value  of  10  Kilo-ohms  is  set  for  the 
tapping  resistors.  The  expression  of  given  in  Appendix  F 
is  used  to  calculate  the  tapping  resistances.  These  values 
are  listed  in  Table  VI . I  in  the  third  column  and  are  used  to 
implement  the  bandpass  filter  using  the  11  taps  of  Reticon 
TAD- 12. 

The  impulse  response  of  this  filter  is  illustrated 
in  Figure  6-3  at  three  different  sampling  frequencies: 

fg  =  50  KHz,  500  KHz,  1.25  MHz 
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TABLE  VI. II.  Theoretical  §  Measured  Coefficients 


k 

A-^  Theor. 

f  =50  KHz 
s 

A]q  Measured 

f  =500  KHz 
s 

f  =1.25  MH 
s 

-5 

0.09863 

0.34444 

0.30508 

0.34483 

-4 

0.07378 

0.33333 

0.209339 

0.06896 

-3 

-0.38389 

-0.5 

-0.66102 

-0.86306 

-2 

-0.59148 

-0.6 

-0.66102 

-0.51724 

-1 

0.28647 

0.48888 

0.49152 

0. 72413 

0 

1. 

1. 

1. 

0.86206 

1 

0.28647 

0.44444 

0.54237 

0.4 

2 

-0.59148 

-0.65555 

-0.76271 

-1. 

3 

-0.38389 

-0.48888 

-0.71186 

-0.34482 

4 

0.07378 

0.31111 

0.13559 

0.51724 

5 

0.09863 

0 . 26666 

0.13559 

0.24138 

The  input 

pulses  used 

are  listed  in  Table 

IV. II  for 

different  sampling 

frequencies 

.  The  tapping  coefficients 

measured  with  the  sample  and  hold  circuit  are  shown  in  Table 
VI.  II.  The  theoretical  and  measured  are  quite  different 
for  some  taps  due  to  the  fixed  pattern  noise  and  the  non¬ 
uniformity  of  the  tap  outputs.  The  maximum  noise  observed 
occurs  between  tap  #3  and  #4  outputs.  Its  peak  amplitude  is 
about  30%  of  A^. 

The  theoretical  magnitude  frequency  response  and 
measured  response  curves  for  three  different  sampling  fre¬ 
quencies  are  plotted  in  Figure  6-4.  The  main  lobes  of  measured 
response  curves  are  rather  similar  to  the  theoretical  one. 
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But  there  are  some  differences  at  the  side  lobes  which  have 
their  peak  magnitudes  as  follows: 

-  15%  of  maximum  amplitude  for  f^  =  50  KHz 

-  35%  of  maximum  amplitude  for  f^  =  500  KHz 

-  51%  of  maximum  amplitude  for  f^  =  1.25  MHz. 

There  are  also  some  wiggles  for  input  frequencies  less  than 

.15  f^. 
s 

2 .  Tapping  Resistors  Adjusted  to  Yield  Better 

Impulse  Response,  Rmin  =  10  Kilo-ohms 

Since  for  k  greater  than  3,  the  coefficients  have 
negligible  values  in  comparison  with  the  noise,  the  adjustment 
is  done  only  on  7  tap  outputs.  Obtaining  the  exact  values  of 
the  theoretical  design  is  very  delicate  because  of  the  inter¬ 
action  between  tap  outputs.  The  adjusted  coefficients  and 
the  corresponding  tap  resistances  are  listed  in  Table  VI. III. 
Some  modification  has  to  be  done  on  the  values  of  some  tap 
resistors  in  order  to  get  better  impulse  response  for  differ¬ 
ent  sampling  frequencies. 

The  amplitude  and  width  of  input  pulses  are  listed 
in  Table  IV. II. 

The  impulse  response  of  this  filter  is  illustrated 
in  Figure  6-5  at  three  different  typical  sampling  frequencies. 
The  maximum  noise  observed  has  its  peak  amplitude  about  19% 
of  Aq. 

The  theoretical  and  measured  curves  of  bandpass  filter 
magnitude  frequency  response  are  shown  in  Figure  6-6.  The 
wiggling  of  the  frequency  response  at  input  frequencies  less 


123 


than  .15  is  much  reduced  for  =  50  KHz,  but  the  side 
lobe  is  still  high  at  15%  of  the  maximum  amplitude.  For  = 
1.25  MHz  and  =  500  KHz,  the  peak  amplitudes  of  side  lobes 
are  increased  to  29%  of  maximum  amplitude. 

TABLE  VI.  III.  Adjusted  Coefficients  ^  Resistors  (Rmin"  10 


k  A^  Rk  (Kilo-ohms) 

f^=50KHz  f^=500KHz  f^=1.25MHz  f  =50KHZ  £  =500KHz  £^=1.25MHz 

s  s  s  s  s  s 


-5 

.07246 

.01587 

.07317 

147. 

147. 

147. 

-4 

.08695 

.04762 

.04878 

198. 

198. 

198. 

-3 

-.20289 

- . 55555 

-.34146 

48.07 

40. 

60.11 

-2 

-.62318 

- .62492 

- . 53658 

23. 

20. 

16. 

-1 

.28985 

.30156 

.26829 

47. 

75. 

210. 

0 

1. 

1. 

1. 

8. 

8. 

6. 

1 

.28985 

.30156 

.26829 

47. 

75. 

9.8 

2 

- . 62318 

- .63492 

- . 53658 

23. 

24.8 

8.5 

3 

-.20289 

- . 55555 

-.34146 

48.07 

40. 

23. 

4 

. 08695 

.04762 

. 07317 

198. 

198. 

198. 

5 

.07246 

.01587 

- .04878 

147. 

147. 

147. 

3 .  Tapping  Resistors  According  to  Theoretical 

Design.  Rmin  =  50  ohms 

A  minimum  value  of  50  ohms  is  set  for  the  tapping 
resistance.  Using  the  formula  described  in  Appendix  F,  the 
tapping  resistances  are  calculated  and  have  the  values  listed 
in  Table  VI . I ,  fourth  column. 

The  pulse  inputs  used  have  the  amplitude  and  width 
listed  in  Table  IV. II.  Figure  6-7  illustrates  the  impulse 
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the  impulse  response  of  the  filter  at  sampling  frequencies: 


f  =  50  KHz,  500  KHz,  1.25  MHZ 
s 


measured 

coefficients  Aj^ 

are  listed 

in  Table  VI. IV. 

TABLE  VI 

.IV.  Measured 

Coefficients 

(R  50  ohms) 

min 

k 

fg=50  KHz 

■^k 

f  =500  KHz 
s 

f  =1.25  MHz 
s 

-5 

.03211 

.03137 

. 82810 

-4 

.05504 

. 07058 

.82062 

-3 

- .37155 

- . 27451 

.  76083 

-2 

- . 58256 

-.40392 

.  74439 

-1 

. 24312 

.22745 

.84603 

0 

1. 

1. 

1. 

1 

. 24312 

.24314 

. 84603 

2 

-.58256 

-.34509 

.73243 

3 

- .37155 

-.38823 

. 73094 

4 

.05504 

.04313 

.83408 

5 

.05504 

.09803 

.81464 

The  maximum  peak  amplitude  of  noise  measured  is  about  18%  of  Aq. 

The  theoretical  and  measured  curves  of  bandpass  filter 
magnitude  frequency  response  are  shown  in  Figure  6-8.  The 
wiggling  effect  of  the  response  curves  are  strongly  diminished 
in  comparison  with  previous  curves  in  Figures  6-4  8  6-6.  The 
amplitudes  of  side  lobes  are  also  reduced  and  have  the  follow¬ 
ing  values : 
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50  KHz 


7.5%  of  maximum  amplitude  for  f^  = 

-  11%  of  maximum  amplitudes  for  f^  =  500  KHz 

-  12.5%  of  maximum  amplitude  for  f^  =  1.25  MHz. 

4 .  Tapping  Resistor  Adjusted  to  Yield  Better 

Impulse  Response.  Rmin  =  50  ohms 

a.  Driving  Circuit  Using  RCA/CMOS  Flip-Flop 

A  minimum  value  of  50  ohms  is  set  for  tap  #6. 

The  tap  outputs  are  adjusted  by  sample  and  hold  circuit  to 
yield  better  impulse  response.  The  tapping  coefficients  and 
resistances  adjusted  are  listed  in  Table  VI. V. 

TABLE  VI. V.  Adjusted  Tapping  Coefficients  §  Resistors 


k  Ak  Rv  (Kilo-ohms) 

f  =50KHz  f  =500KHz  f  =1.25  MHz  f  =50KHz  f  =500KHz  f  =1.25  MHz 
s  s  s  s  s  s 


-5 

.16733 

. 08532 

.01234 

46.2 

46.2 

46.2 

-4 

. 06374 

.07508 

. 03086 

63.45 

63.45 

63.45 

-3 

- .37450 

- .33673 

-.25308 

6.25 

4.05 

4.05 

-2 

- .60557 

-.59183 

-.59876 

1.2 

1.2 

.  5 

-1 

. 28287 

.23849 

.22222 

23.023 

25.023 

10. 

0 

1. 

1. 

1. 

.  05 

.05 

.  05 

1 

. 28287 

.23849 

.23457 

6.1 

6.1 

4.1 

2 

-.60557 

-.59183 

-.59876 

1.083 

1.083 

1.08 

3 

-.37450 

-.33673 

- .23457 

4.  083 

4.55 

1.083 

4 

.06374 

.07508 

. 04321 

63.45 

63.45 

63.45 

5 

.16733 

.08532 

.02469 

46.2 

46.2 

46.2 

The 

input  pulse  used 

have  the 

amplitudes  and 

widths 

listed  in 

Table  IV. II.  Figure  6-9  illustrates  the  impulse  response  of 
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the  filter  at  three  different  typical  sampling  frequencies. 
The  maximum  noise  has  its  peak  amplitude  about  15%  of  A^. 

The  theoretical  and  measured  curves  of  bandpass 
filter  magnitude  frequency  response  are  shown  in  Figure  6-10. 
The  wiggling  does  not  appear  at  the  main  lobe  but  the  ampli¬ 
tudes  of  the  side  lobes  are  s.till  high: 


-  15% 

of  maximum 

amplitude 

for 

£3  =  500 

KHz 

-  10% 

of  maximum 

amplitude 

for 

£3  =  500 

KHz  §  1.25  MHz 

b.  Driving  Circuit  Using  INSELEX  CMOS/SOS  Flip-Flop 

A  new  circuit  with  faster  flip-flop  is  made  and  a 

minimum  value  of  50  ohms  is  set  for  tapping  resistance.  The 

adjustment  is  done  for  f^  equal  to  500  KHz;  the  coefficients 

and  resistances  obtained  are  used  also  for  f  of  50  KHz  and 

s 

1.25  MHz.  Table  VI. VI  lists  these  coefficients  and  resistances. 

TABLE  VI. VI  Adjusted  Coefficients  ^  Resistances 
k  A]^  R]^  (Kilo- ohms) 


-5 

. 03536 

46.2 

-4 

.01571 

63.45 

-3 

- .40078 

3.8 

-2 

-.61689 

1.56 

-1 

.29194 

9.023 

0 

1. 

.05 

1 

.29273 

2.32 

2 

- . 61297 

1.263 

3 

-  .40078 

5.000 

4 

.06287 

63.45 

5 

.07269 

46.  2 
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The  pulse  inputs  listed  in  Table  IV. II  are  used.  Figure  6-11 
illustrates  the  impulse  response  of  the  filter  at  sampling 
frequencies . 

f  =  50  KHz,  500  KHz,  1.25  MHz. 
s  *  * 

The  maximum  peak  amplitude  noise  observed  is  about  21.6%  of 
^o* 

The  theoretical  and  measured  curves  of  bandpass 
filter  magnitude  frequency  response  are  shown  in  Figure  6-12. 
There  is  no  wiggling  at  the  main  lobe  region  and  some  improve¬ 
ment  is  obtained  at  the  side  lobes,  which  have  the  following 
amplitudes : 

6%  of  maximum  amplitude  for  f^  =  1.25  MHz 

-  10%  of  maximum  amplitude  for  f^  =  500  KHz 

-  14%  of  maximum  amplitude  for  f^  =  50  KHz. 

Using  the  computer  program  shown  in  Appendix  E, 

the  measured  magnitude  frequency  response  for  f^  equal  to 
500  KHz  and  the  response  corresponding  to  the  measured  impulse 
response  listed  in  Table  VI. VI  are  plotted  in  Figure  6-13  with 
DB  scale  on  Y  axis  and  normalized  frequencies  on  X  axis.  The 
main  lobes  of  the  two  curves  are  the  same  but  there  are  some 
differences  at  the  side  lobe  regions.  This  shows  that  the 
technique  using  the  sample  and  hold  circuit  to  adjust  the 
tapping  resistor  to  get  better  impulse  response  is  quite 
accurate  and  very  useful  to  the  filter  implementation. 
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5 .  Distortion  Measurements  of  Bandpass  Filter 

For  the  case  described  in  V.B.3  where  the  tapping 

resistances  are  set  according  to  theoretical  design  with 

R  .  equal  to  50  ohms,  the  second  harmonic  distortion  de- 
min  ^  ’ 

fined  by: 


DB 


20  1  P  2nd  harmonic  voltage 
fundamental  voltage 


are  illustrated  respectively  in  Figures  6-14  and  6-15,  for 
sampling  frequency  of  500  KHz  and  three  different  peak  to 
peak  input  voltages: 


V.  =  3V,  0.8V,  O.IV 
in  *  ’ 

From  the  Figure  6-14  where  the  input  frequency  is  equal  to 
50  KHz,  the  second  harmonic  and  clock  noise  at  85  KHz  have 
the  relative  values  expressed  in  percentage  of  fundamental 
as  shown  in  Table  VI. VI I. 


TABLE  VI. VI I  Second  Harmonic  8  Clock  Noise 


V.  ppk 
in 

3  Volts 

. 8  Volts 

.1  Volt 

2nd  Harmonic 

16.7% 

6.5% 

6.  8% 

Clock  Noise 

5.8% 

14.9% 

88.8% 

When  the  amplitude  of  the  input  voltage  decreases,  the  rela¬ 


tive  value  of  second  harmonic  decreases  but  the  relative  value 


of  clock  noise  increases.  In  the  Figure  6-14c,  one  can  observe 


many  spurious  responses  which  occur  at  63  KHz,  85  KHz,  112  KHz, 


126  KHz,  and  are  due  principally  to  the  clock  noise  and  fixed 


pattern  noise. 

f 

} 
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oeo  S2C  020  OTC  CTO  5C0 


FIGURE  6.1  Amplitude  Frequency  Response  of  Bandpass  Filter 
Length  =  11  /  Hamming  Window 

X  -  Scale  =  .2  units/inch  for  normalized  frequency 
Y  -  Scale  =  10  units/inch  for  amplitude  frequency 
response. 
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FIGURE  6*2  Amplitude  Frequency  Response  of  Bandpass  Filter  . 

Len^h  =  11  /  Hamming  Window 

X  -  Scale  »  .2  units/inch  for  normalized  frequency 
Y  —  Scale  »  10  units/inch  in  DB  for  amplitude  frequency 
response. 
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50  KHz  500  KHz  1.25  MHz 
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FIGURE  ^,4  Beindpass  Filter  Magnitude  Frequency  Response, 

Tapping  Resistors  According  to  Design,  R  .  =10  Kilo-ohms 


50  KHz  500  KHz  1.25  MHz 
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FIGURE  6,6  Bandpass  Filter  Magnitude  Frequency  Response, 

Tapping  Resistors  Adjusted  to  Yield  Better  Impulse  Response,  R^^^  -  10  Kilo-ohms 
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FIGURE  6.8  Bandpass  Filter  Magnitude  Frequency  Response. 

Tapping  Resistors  According  to  Design.  R  “  50  ohms 


50  KHz  -500  KHz  1.25  MHz 
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FIGURE  6.10  Bandpass  Filter  Habitude  Frequency  Response. 

Tapping  Resistors  Adjusted  to  Yield  Better  Impulse  Response.  R  -50  ohms 
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.FIGURE  6*12  ijandpass  Filter  Magnitude  Frequency  Response. 

Tapping  Resistors  adjusted  to  Held  Oetter  Impulse  Response.  R^^  =  50 


900-  OTO-  910-  OZO-  9ZC-  OrO-  ^jZO-  Cl-O- 


FIGURE  6.13 


Bandpass  Filter  Magnitude  Frequency  Response. 

Tapping  Resistors  Adjusted  Using  Impulse  Response. 

Rmin  =  50  ohms.  Q  =  ,250,.  J 
X^*^-^ScaIe  =  .1  units/inch  £c 

Y  -  Scale  =  5  units/inch  in  DB 
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FIGURE  6.15  Second  Harmonic  Distortion  of  Bandpass  Filter  with  Tapping  Resistors  According 

to  Design  and  50  ohms  for  Different  Input  Signal  Amplitudes,  fga^^p  -  500  KHz 


6.  Observations 


The  characteristic  feature  of  the  filter  realized 
by  the  electronic  tapped  analog  delay  line  is  that  one  can 
scan  the  center  frequency  of  the  bandpass  filter  to  the  left 
or  right  by  simply  decreasing  or  increasing  the  clock  fre¬ 
quency.  Figure  6-16  illustrates  this  feature  for  Q  =  constant 
at  two  different  sampling  frequencies: 


FIGURE  6-16.  Bandpass  Filters  Amplitude  Frequency 
response  at  2  different  sampling 
frequencies  f  =500  KHz  §  f  =835  KHz. 
Scale :  H  =  .  2^nsec/Div  V  ^ .  1  V/Div. 


Another  observation  is  that  the  transversal  filter 
implementation  is  not  very  sensitive  to  high  specifications 
for  moderate  number  of  taps.  The  tapping  coefficients  are 
almost  the  same  for  different  values  of  Q  required.  The 
impulse  response  of  bandpass  filters  for: 

f^  =  0.2  fg  (f^  =  center  frequency) 

Q  =  50,  250,  500 
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using  11  taps  is  listed  in  Table  VI. VIII,  and  using  23  taps 
in  Table  I. IX. 


TABLE  VI. VI I I  Tapping  Coefficients  (11  Taps) 

k  Ak 

Q=50  Q=250_ Q=500 


0 

1. 

1. 

1. 

1 

0.28644 

0.28647 

0.28640 

2 

-0.59140 

-0.59148 

-0.59148 

3 

-0.38382 

-0.38389 

-0.38395 

4 

0.07375 

0.7378 

0.07378 

5 

0.09857 

0.09863 

0.09863 

TABLE 

VI. IX  Tapping 

Coefficients 

(23  Taps) 

k 

Q=50 

Ak 

Q=250 

Q=500 

0 

1. 

1. 

1. 

1 

0.30373 

0.30377* 

0.30370 

2 

-0.75475 

-0.75486 

-0.75486 

3 

-0.69072 

-0.69086 

-0.69095 

\ 

4 

0.23217 

0.23227 

0.23227 

5 

0.63317 

0.63357 

0.63358 

6 

0.15702 

0.15717 

0.15718 

7 

-0.31184 

-0.31222 

-0.31224 

8 

-0.22189 

-0.22225 

-0.22226 

9 

0.05648 

0.05660 

0.05660 

10 

0.11777 

0.11807 

0.11808 

11 

0.02596 

0.02604 

0.02605 
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VI.  CONCLUSION 


Digital  non- recursive  filters  have  been  extensively 
developed  for  the  implementation  of  "finite  impulse  response" 
(FIR)  filters  using  software  on  digital  computers.  In  recent 
years,  new  CTD  (charge  transfer  device)  tapped  delay  lines 
make  it  possible  to  implement  real  time  FIR  filters  using  these 
hardware  devices.  However,  because  the  signal  is  sampled 
analog  instead  of  digital  and  because  the  device  principle 
and  technology  are  new,  careful  investigation  is  needed  for 
the  successful  development  of  these  hardware  sampled  analog 
signal  processors.  In  this  thesis,  the  new  Reticon  TAD-12 
tapped  delay  line  with  external  programmable  tap  weights  is 
investigated  in  detail.  In  addition  to  a  thorough  evaluation 
of  its  device  performance,  the  properties  of  three  sampled 
analog  FIR  filters  using  TAD- 12  are  investigated.  They  are 
the  prewhitening,  dewhitening  and  bandpass  filters. 

Both  the  frequency  response  of  TAD- 12  devices  have  been 
evaluated  for  sampling  frequency  up  to  2.5  MHz.  It  was  found 
that  the  output  showed  the  following  limitations: 

1.  Frequency  response  limitation: 

Tap  outputs  generally  roll-off  faster  than  predicted 
by  the  effect  of  sample  and  hold.  The  roll-off  depends  on 
the  location  of  the  tap  and  its  tapping  resistance. 

2.  Time  response  limitation: 
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Output  at  one  tap  generally  takes  more  than  the  delay 
period  between  adjacent  taps  to  settle  down  and  affect  the 
output  of  the  following  taps.  The  settling  time  depends  on 
the  location  of  the  tap,  its  tapping  resistance  and  the 
sampling  frequency. 

3.  Non-uniformities: 

Tap  outputs  are  not  uniform  among  the  taps.  The  non¬ 
uniformity  varies  with  the  location  of  the  tap,  its  tapping 
resistance  and  signal  frequency. 

4.  Loading  effect: 

Tap  outputs  vary  with  tapping  resistance.  However, 
it  should  be  pointed  out  that  these  limitations  are  large  at 
higher  sampling  frequencies  and  signal  frequencies.  At 
sampling  frequencies  less  than  100  KHz  and  signal  frequencies 
less  than  10  KHz,  TAD-12  behaves  satisfactorily. 

The  properties  of  three  sampled  analog  FIR  filters  using 
TAD-12  have  been  evaluated  following  three  design  procedures. 

1.  Tap  resistances  are  selected  using  the  same  design 
procedure  developed  for  digital  non-recursive  filters. 

2.  Tap  resistances  are  corrected  using  the  frequency  response. 

3.  Tap  resistances  are  corrected  using  the  impulse  time 
response . 

It  was  found  that  the  device  limitation  causes  considerable 
deviation  of  the  filter  performance  from  theoretical  prediction 
especially  when  the  sampling  frequency  is  high.  Furthermore, 
the  deviation  is  relatively  smaller  in  bandpass  filter  and 
becomes  worse  in  dewhitening  filter  followed  by  prewhitening 
filter.  This  trend  is  the  result  of  the  relative  magnitudes 
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of  the  tap  outputs  with  respect  to  the  noises  caused  by  the 
device  limitations. 

These  problems  encountered  in  the  Reticon  TAD-12  tapped 
delay  line  just  showed  that  it  is  still  in  its  developing 
stage.  Improved  results  can  be  expected  when  the  causes  of 
the  device  limitations  are  better  understood  and  corrected. 
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APPENDIX 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 


c 

c 


4:^:  ^  ^  :4c:^:{c;jc:{c:4c3)(:!c  :4i :4c :41:41:4c 

FFOGFAP  OF  C6SIGNING  NQNRECURSIVE  PREWHITENING  FILTER  BY 
VINCCW  METHOD 

« :4t « :(i  3jc  :4.  :4c  :4i  :4!  cjt  34c  :jt  :4:  *  :4c  34c  :(t  :4c  !4t  4=  54!  ♦♦*♦♦*  * 


THIS  PROGRAM  CALCULATES  AND  PLOTS  AMPLITUDE  ANC  PHASE  OF 
FIR  CIGITAL  FILTER  WITH  TRANSFER  FUNCTION 
H(Z)  =  C(0)+C(1  )*Z*=!'(-l)  +  C(2)*Z=<'*(-2)+....+C(M  )*2=i'*(-M) 
C(M)  FINITE  IMPULSE  RESPONSE 

A  SLOPE  CF  FREQUENCY  AMPLITUDE  RESPONSE  FRCM  FO  TO 

HIGHEST  FREQUENCY  FP 

R  CONSTANT  FREQUENCY  AMPLITUDE  RSEPONSE  FRCM  0  TO  FG 

T  =  1.  (NORMALIZED  DELAY  UNIT)  AND  FP=0.5=<'(  1/T  )  AND 
£  =  FQ/FP 


COMPLEX  ZINV,G,P 
COMPLEX  CEXP.CMPLX 

DIMENSION  F(150),C(30) ,W(30  ),0(30)tDD(30),CW(20), 
lAMPdSC  )tPHASE(  150) 

INTEGER’»4  ITB  ( 12  )/ 12*0/ 

REAL*A  RTB(28)/28*0.0/ 

CATA  R/0.1/,A/1./ ,E/119.048/,N/5/ 

PI=2. 14159265 

CCNVRT=180./PI 

M=2*N+1 

WFITE(6,500)  M 

5QC  FORMAT! ix, 'LENGTH  OF  FILTER  IS:', 14///) 

TO  COMPUTE  FINITE  IMPULSE  RESPONSE  ANC  WINDOW  FUNCTION 


W  (N*l)=l. 

0 (N+1 )=R+A/(2*PI )*( .5-l./(2.*E) )**2 
CO  600  K=1.N 

C  (K)=A*  (CCS(K*Pn-CCS  (K*PI/E)  )/  (  2.*PI**3*K  = 


C  (N  + 


C  (N  +  l+K  )  =  C(N+1-K) 


DC  (K 


-K )=D(K) 


*2  ) 


=0.5^+0.46*CQS(2.*PI*K/M) 


W(N  +  1-K  )=0D(K) 

W (N  +  l+K  )  =  W(N+1-K) 

600  CONTINUE 

WRITE (6,650) (C( I) ,I=1,M) 

650  FCFMAT(//2X, 'COEFFICIENTS  A  RE : ' /8 ( IX , FIO. 5 )  ) 
CMAX=A6S(C(N+1 ) ) 

DC  700  1=1, M 

IF(ABS(C(I )).GT.CMAX)  CMAX=A8S(C(I)) 

700  CONTINUE 

CO  750  1  =  1, M 

NORMALIZE  FIR 

C  (I)  =  C(  D/CMAX 
750  CONTINUE 

WFITE(6,800)(C( I) ,I=1,M) 

800  FCPMAT( //2X, 'NORMALIZED  COEFFICIENTS  AR E : • / £ ( IX, F 10. 5 ) 1 
WRITE  (6,900)  (W(I),I=1,M) 

900  FCFMAT(//2X, 'HAMMING  WINDOW  ARE  :  « /8( IX, FIO. 5)  ) 

DC  9  1=1, M 
9  CW(I)=C (I )*W(I ) 

WFITE(6,1000) (CW( I  )  ,I  =  1,M  ) 

1000  FORMAT! //2X,' COEFFICIENTS  IMPROVED  BY  HAMMING  ARE' 
l/8( 1X,F10.5) ) 

WFITE(6,50) 

50  FORMAT! ///5X,' FREQUENCY ',8X,* AMPLITUDE' , 12X  ,  • PHAS 5  ' ///  ) 
CO  10  J  =  1,1Q1 


F  IS  FREQUENCY  IN  FRACTIONS  OF  NYQUIST  FREQUENCY 
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c 

FU)=0,C05«PL0AT(  J-1) 
2INV=CEXP(CMPLX(0.,-2.*FC J)*PI  )  ) 

P  =  C^.PLX{1.,0.) 

C-sCf'PLXCO.tO.) 

CC  6  1  =  1, M 
G=G4p*CW(I ) 

6  P=F*Z!NV 

Af'P  (J)  =  CABS(G) 

X=PEALCG) 

V=AIMAG(G» 

IF(ABS(X).GT.l.E-8)  GO  TO  60 
FF4SE(J)=SIGN{90.,Y> 

GC  TO  8 

60  PFASE(J)=ATAN2(Y,X)*C0NVRT 
8  UPITE16,20)  F( J) , ANP(J» ,PHASE( J  ) 

20  FCPMAT(5X,F6.3,4X,E16.8,4X,E16.e> 

10  CCMINUE 

ECLIVALENCE(TITLE,RTB{5) ) 

RE^L^e  TITLE(12) 

PEAC(5,300)  TITLE 
300  FCFMAT(6A8) 

C^LL  CRAWP(101,F,AMP,IT8,RTe) 

PE4C(5,200)  TITLE 

C^LL  DRAWPdOl, F, PHASE, ITB,  RTB) 

STOP 

ENC 

//GC.SYSIN  CC  * 

AMPLITUCE  FRECUENCY  RESPONSE  OF  PREWHITENING 
FILTER  /  LENGTH  =  11  /  HARMING  WINDOW 
PHASE  RESPONSE  OF  PREWHITENING  FILTER 
LENGTH  *  11  /  HAMMING  WINDOW 
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APPENDIX  E 


C  FFGGFAM  OF  CESIGNING  DWF  FROM  GIVEN  IMPULSE  RESPONSE  OF 
C  FWF  EY  FRECUENCY  SAMPLING  METHOD 

C  TFE  FPESCRIBEC  CHARACTERISTIC  AND  DESIGNED  FRECUENCY 
C  CFAP/CTEPISTIC  OF  DWF  ARE  PLOTTED  ON  THE  SAME  GRAPH 

C  4=  ’S' ^  4:  =!'=!' ’S' T  =4=  =4=  * ’i' * ’i' # 


THIS  PROGRAM  CALCULATES  AND  PLOTS  AMPLITUDE  ANC  PHASE  OF 
FIR  DIGITAL  FILTER  WITH  TRANSFER  FUNCTION 
H(2»  =  C(0  )+C(l  )*Z**(-1  )  +  C(2 )*Z**(-2)+....+C(M  )«Z**( -M ) 
C(M)  FINITE  IMPULSE  RESPONSE 


CCHFLEX  ZINV,G,P»H(150) 

COMPLEX  CEXP.CMPLX 

DIMENSION  F(150),C(30),W(30),D(20),0D(30)»CW(30), 
lAMP (1 50 )» PHASE ( 150 )»HH( 150) 

INTEGER=»4  ITB(  12  )  /  12*0/ 

REAL*4  RTB(28  )/28*0.0/ 

FACT0R=1.12 
MM=0 
N  =  5 

PI=3.1415S265 
CCNVRT= 180-/PI 
M=2*N+1 

WFITE<6,500)  M 

500  FCFMAT( IX, 'LENGTH  OF  FILTER  IS:',IA///) 

READ(5, 1) (C(I) ,I=1,M) 

1  FCRMAT( EF1C.5) 

610  WRITE (6,650) (C( I ) ,1=1, M ) 

650  FCFMAT( //2X, 'COEFFICIENTS  ARE : ' /8 ( IX , FIO .5  )  ) 

CCMPLTE  HAMMING  WINDOW 

W  (N+l)=l. 

DC  6C0  l<  =  l,N 

CC(K)=0.54+0.46*C0S (2.*PI*K/M) 

W(N  +  1-K  )  =  0D(K ) 

W  (N  +  l+K )=W (N+l-K ) 

60C  CONTINUE 

CMAX=ABS(C(N+1) ) 

DC  700  1  =  1, M 

IF(ABS( C(I ) ).GT.CMAX)  CMAX= ABS ( C ( I ) ) 

700  CONTINUE 

DC  750  1  =  1  ,M 

NORMALIZE  FIR 

C (I )=C( i)/CMAX 
750  CONTINUE 

WRITE (6,800)(C( I ) ,I=1,M) 

800  FCRMAT( //2X, 'NORMALIZED  COEFFICIENTS  ARE : * /8( IX, F10.5)  J 
DC  5  1=1, M 

<  cw(r)=c (I )*w(i ) 

WPITE(6,1000)(CW( I )  ,I  =  1,M) 

1000  FCRMAT( //2X,'CCErFICIENtS  IMPROVED  BY  HAMMING  ARE' 
l/e(lX,F10.5) ) 

WFITE(6,50) 

50  FC F MAT (///5X,'FREQUENCY',8X,' AMPLITUDE' , 13X , • PHAS E ' / //  ) 
CO  10  J=l,101 

F  IS  FREQUENCY  IN  FRACTIONS  OF  NYQUIST  FREQUENCY 

F(  J  )=0.Q05*FL0AT(  J-1) 

ZINV=CEXP(CMPLX(0.,-2.*F( J )*PI ) ) 

F=CMPLX(1.,0.) 

G=CMFLX (0.,0.) 
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or>o  ooooo  ooo  ooo 


11 

e 


cc  6  1=1, M 

IF(MM.£C.l )  GO  TO  11 
G=G+P*CW(I ) 

GC  TO  6 

G  =  G4P*CW(n*FACT0R 

F=F*ZINV 

^^'F  (J)=CABS(G) 


FPECteNCY  SAMPLE  POINTS  OF  DWF 

l-(J)  =  l./G 
HK  J)=C4BS(H(J)  ) 

X=REAL(G)- 

Y=AIMAG{G) 

IF(AeS(>().GT.l.E-8)  GO  TO  60 
PF^SE(J)=SIGN(90.,Y) 

GC  TO  8 

6C  PF^SEt J)=ATAN2(Y,X)*C0NVRT 
8  WPITE<6,20)  F( J) ,AMP( J) ,PHAS£( J) 
2C  FCFMAT(5X,F6.3,4X,£16.8,4X, E16.8) 
10  CONTINUE 

ECU IVALENCE (TITLE, RTB(5) ) 

RE4L*8  TITLE(12) 

IF(MM.EC.l)  GO  TO  301 
GC  TO  302 
301  ITE(1)=2 
ITE(2)=0 

C4LL  CRAWPdOl,  F,AMP,  ITB,RTB) 

202  L=N+1 


CF  CWF 


1  =  1, L 


HOC 


1500 


300 


CC  1500 
l<  =  I-l 

6G»CMPLX(0.,0. ) 

CC  1100  J=l,101 

2=CEXP(CMPLX(0.,2.*PI*F(J )#K) ) 
GG=6G+2#H( J) 

CCNTINUE 
C  (n  =  66/101. 

C (M  +  l-I )  =  C(I  ) 

CCNTINUE 


NN=MM+1 
IF(MM.EC.2» 
PE4C{5, 200) 
FCFMAT(6A8) 
ITE(1)=1 
ITE(2)=1 


GO  TO 
TITLE 


2100 


FLCT  THE  PRESCRIBED  FREQUENCY  CHARACTERISTIC  CF  CWF 

LSING  INVERSE  DISCRETE  FCURISR  TRANSFORM  TO  CALCULATE  THE 
IMPULSE  RESPONSE  OF  CWF 

CALL  CRAWP(101,F,HH,ITB,RTB ) 

GC  BACK  TO  PLOT  THE  AMPLITUDE  FREQUENCY  RESPONSE  OF  DWF 

610 


GC  TO 
2100  STCF 
ENC 

//GC.SYSIN  CD 

-0.00091  0.0000 

-0.21^27  C.O 
AMPLITUDE  FRECUENCY 
FILTER  /  LENGTH  =  11 


-0.01219 
-0.01219 
RESPONSE  OF 
/  HAMMING 


0.0000 
0.0 
DEWHITEMNG 
WINDOW 


-0.21A27 

-0.00091 


1.0000 
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APPENDIX  C 


Q  #  :»5!!4e  :je :{:  sis  *  4:  sit Sr  tS  *  nS ’Jt  3*;  ^  5>  ♦  Jf! ^  4  ^  ^  * 

C*  FFCGFAM  OF  DESIGNING  NONFECURSIVE  BANDPASS  FILTER 

«  4t  J»  sji  *  ♦  si  ♦ :» :|t  » s{s  «  «  ♦  s{e  Jfe :{:  jjt «  ♦  sjc  sit  ♦  ♦  ♦  Sj:  3}:  ♦  Jjt  #  #  3)1  #  sjs  sis  *  *  ♦  *  s}t « sic  *  *  si.  ♦  ♦  «  sit « s{.  ♦  3*  3*  sji  sit 


TFIS  PROGRAM  CALCULATES  AND  PLOTS  AMPLITUDE  AND  PHASE  OF 

FIR  DIGITAL  FILTER  WITH  TRANSFER  FUNCTION 

K2)  =  C(0)+C(l)*Z**(-l)  +  C(2)=!‘Z**(-2)+....+C(M  )«Z**(-M) 

C(M)  FINITE  IMPULSE  RESPONSE 

A  AMPLITUDE  OF  BANDPASS 

FI  e  F2  I  CRITICAL  FREQUENCIES 


CCMFLEX  ZINVfGfP 
CCMFLEX  CEXPrCMPLX 

DIMENSION  F(150)  tCOO)  ,W(  30),D<  20)  ,DD(30)  fCVtOC)  , 
lAMPdSO  ),PHASE(150) 

INTEGER*4  ITB<  12  ) /12’i'O/ 

REAL*4  RTB(28)/28*0.0/ 

DATA  A/l./,Fl/0.1980/,F2/0.2020/,N/5/ 

PI=2. 14159265 

CCNVRT=180./PI 

M=2*N+1 

WP1TE(6,500)  M 

5QC  FCRMATdX, 'LENGTH  OF  FILTER  IS:*, 14///) 


W  (N+l)  =  l. 

C  (N  +  1)=2.*A«(F2-F1) 

DC  6C0  K=1,N 

D  (K  )  =  A*  (SIN(2.’'‘PI*K’!'F2)-SIN  (2.*PI*K*F1)  )/(K’»PI  ) 

C  (N  +  l-K  )=D(K) 

C  (N+1+K)  =  C(N+1-K) 

DC  (K)  =0.54+0. 46*C0S  (2.*PI=>K/M) 

W(N  +  1-K  )=CC(K) 

W (N  +  l+K  )  =  W(N  +  1-K) 

600  CONTINUE 

WRITE  (6,650) (C( I ) ,I  =  1,M) 

650  FCRMAT(//2X, 'COEFFICIENTS  ARE : ' / 8 ( IX, F 10. 5 )  ) 
CMAX=ABS(C(N+1) ) 

DC  TOO  1=1, M 

IF(ABS(C (I  )  ).GT.CMAX)  CMAX=ABS ( C ( I  ) ) 

700  CONTINUE 

DC  750  1=1, M 

NORMALIZE  FIR 

C  (I  )  =  C(I )/CMAX 
75C  CONTINUE 

WRITE (6,800) (C( I) ,I=1,M) 

8QC  FCRMAT(//'2X, 'NORMALIZED  COEFFICIENTS  AR  E :  '  /  E  (  IX,  F  10 . 5 )  ) 
WRITE  (6,900)  (W(I),I=1,M) 

900  FORMAT (//2X, 'HAMMING  WINDOW  ARE :  ' /8 ( IX , F 10. 5 )  ) 

DC  9  1=1, M 
9  CW(I)=C(I  )*W(I ) 

WPITE(6,1000)(CW( I ) ,I=1,M) 

1000  FCFMAT(//2X, 'COEFFICIENTS  IMPROVED  BY  HAMMING  ARE* 
l/e(lX,F10.5)) 

WRITE(6,50) 

50  FORMAT ( ///5X , • F RE QUENCY * , 8X , • AMPL ITUDE ' , 13X , • PHAS E * ///  1 
DC  10  J  =  l,  101 

F  IS  FREQUENCY  IN  FRACTIONS  OF  NYQUIST  FREQUENCY 

F(J)=0.01=»FL0AT(J-1) 

ZINV=CE)(P(CMPLX(0.,-2.=*‘F(  J)*PI  )  ) 

F=CMPLX(1.,0.) 

G=CMPLX (0.,0.) 
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CC  6  1=1, M 
G=G+P*CW(I ) 

F  =  F=»ZINV 
^NP(J)=CABS(G) 

X=FEAL(G) 

Y=AIMAG<G) 

IF(ABS(X).GT.l.E-8)  GO  TO  60 
PI-/SE{J)  =  SIGN(90.,Y) 

GC  TO  8 

PFASE( J )=ATAN2( Y,X)«CONVRT 
V»PITE(6,20)  F(J),A^'P(J),PHASE(J) 
FOPFAK  5X,F6.3,4X,E16.8,4X,E16.8 ) 
CCNTINUE 

ECUIVALENCE(TITLE,RTB(5) ) 

REAL*8  TITLE(12) 

REAC(5,20C)  TITLE 
FCF^AT(6A8) 

CALL  DRAWP(101,F,AMP,ITB,RTE) 
PEAC(5,30C)  TITLE 
CALL  DRAWP(10l,F, PHASE, ITBfPTB) 
STOP 
E^C 

//GC.SYSIN  DC  * 

AMPLITUDE  FREQUENCY  RESPONSE  OF  BANDPASS 
FILTER  /  LENGTH  =  11  /  HAMMING  WINDOW 
PHASE  RESFCNSE  OF  BANDPASS  FILTER 
LENGTH  =  11  /  HAMMING  WINDOW 


6C 

E 

20 

10 


iOQ 
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APPENDIX  D 


C  THIS  PROGRAM  CALCULATES  IMPULSE  RESPONSE  GIVEN  AMPLITLDE 
C  FPEQLENCY  RESPONSE  OF  FILTER 
COMPLEX  CEXP.CMPLX 

CIMENSICN  F  (150)  ,  COO),  WOO  ),CW  (30 ,00(301  ,F(  150) 

N  =  5 

FI=3.1A15S265 

F=2*N+1 

WFITE(6,5C0)  M 

5CC  FCPMAT( IX, 'LENGTH  OF  FILTER  IS: *,I4//) 

COMPUTE  HAMMING  WINDOW 

W(N-H)=1. 

DC  EOO  K=1,N 

CC  (K)=0.5  4+0.46«C0S  (2.  =<'PI  *K/M  ) 

W(N+1-K)=00(K) 

W(N+1+K  )  =  W(N-H-K) 

600  CONTINUE 

PEAC(5,'?00)  (H(  J)  ,J  =  1,101) 

700  FGFMAT(6F10*5) 

WRITE(6,750)  (  H(  J  )  ,  J*  1,  10  1 ) 

75C  F0RMAT(//2X,'DATA  ARE  :',///( 8 ( 1 X  ,  F10.5 ),/) ) 

NCFMALIZE  KJ) 

FMAX=H( 69) 

CC  2000  J=l,101 

IF(H( J) .GT.HMAX)  HMAX=H(J) 

2000  CONTINUE 

CC  2100  J=l,101 
H(J)=H( Jl/HMAX 
2100  CONTINUE 

LSIN6  INVERSE  DISCRETE  FOURIER  TRANSFORM 
TC  CALCULATE  THE  IMPULSE  RESPONSE 


L=M 

CC  1500  1=1, L 
K  =  I-1 

GC-=CMPLX(0,,0.  ) 

DC  1100  J  =  l,101 
F(J)=0.01*FLOAT(J-1) 

Z=CEXP(CMPLX(0.,2.*PI«F(J)*K) ) 

GG  =  GG  +  Z=»H(J) 

HOC  CONTINUE 

C  (I  )  =  GG/1C1. 

1500  CONTINUE 

WPITE(6,650)(C(I) ,I=1,M) 

650  FCRMAT(//2X, 'COEFFICIENTS  ARE : ' /8 ( IX , F 10. 5 )  ) 

NORMALIZE  FIR- 

CMAX=ABS(C(N+1 ) ) 

CC  800  1=1, M 

IF(ABS( C(I n.GT.CMAX)  CMAX=ABS ( C ( I ) ) 

800  CONTINUE 

CC  900  I=1,M 
C(I  )  =  C(  D/CMAX 
900  CONTINUE 

WRITE (6,850) (C( I )  ,1  =  1 ,M) 

850  FCFMAT(//2X, 'NORMALIZED  COEFFICIENTS  AR E  :  •  / £ ( IX , F 10. 5 ) 
CC  9  1=1, M 
9  CW(I)=C(I)=f=W(I) 

WRITE (6, 1000 ) (CW( I )  ,I  =  1,M  ) 

1000  FCRMAT( //2X, 'COEFFICIENTS  IMPRCVEC  BY  HAMMING' 
l/8(  1X,F10.5)  ) 

S7CF 

ENC 
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APPENDIX  E 


C  PPCGPAM  FOR  PLOTTING  100  CALCULATED  AND  MEASURED  POINTS 
C  CF  ANY  NONPECURSIVE  FILTER  WITH  GIVEN  IMPULSE  RESPONSE 

C  :!i  ♦  41  4  *  :Jr  jjt  :ji  :»  4t  ♦  :^c  sjc  jjc  :jt:Jc  *  sjt  *  *  *  *  5jt  ^  :jc  :jt  :jc :{;  jjt  *  4e  :je :» jjt  ;|t «  4e  *  *  *  ♦  *  :{c :(!  ♦ :(!  jjc 


This  PROGRAM  CALCULATES  AND  PLOTS  AMPLITUDE  AND  PHASE  QF 
FIR  DIGITAL  FILTER  WITH  TRANSFER  FUNCTION 
K2)  =  C{0  )+C(l  )*Z**(-l)  +  C(2)*Z**(-2)+....+C(M  )*Z*=«'(-M) 
C(M)  FINITE  IMPULSE  RESPONSE 


CCNPLEX  ZINV,G,P 
CCNFLEX  CEXP.CMPLX 

DIN  ENSIGN  F(150) ,C(30) ,AMP( 150) , PHASE (150) ,8(150) , 

16  6(150  )  ,B A (150), AMPL( 150), AMPLT( 150) 

1AVPL( 150) ,AMPLT(150) 

INTEGER^4  ITB ( 12  ) /12*0/ 

REAL*4  RTB(28  )/28*0.0/ 

PI=2. 1^159265 
N  =  ll 

CCNVRT=180./PI 
PEAC(5, 500)(C( I) ,I=1,M) 

50C  FCPNAT(8F10.5) 

REAC(5,500)  (B(K) ,K=1,101) 

NPITE(6,850)  (B(K) ,K=1, 101) 

85C  FCPMAT(//2X,' DATA  ARE : • ,/// (8 ( IX  ,  F10.5 ) , / )  ) 

BMAX=B( E9) 

CC  1000  J=l,101 

1000  IF(e( J).GT.BMAX)  BMAX=B(J) 

CC  1100  J=l,101 
EE(J)=e(J)/BMAX 
1100  eA( J)=20.*ALOG10(B6(J ) ) 

WRITE (6,80C)(C( I),I=1,M) 

800  FCFMAT(//2X, 'NORMALIZED  COEFFICIENTS  AR E :  •  / €( IX, F 10. 5 )  ) 
WRITE(6,50) 

50  FC P MAT (///5X, 'FREQUENCY', 8X , • AMPL ITUDE ’ , 13X ,' PHASE • /// ) 
CC  10  J=l,101 

F  IS  FREQUENCY  IN  FRACTIONS  OF  NYQUIST  FREQUENCY 

F( J)=0.005*FLOAT( J-1) 

ZINV=CEXP(CMPLX(0.,-2.=!‘F(J  )*PI  )  ) 

F=CNPLX(1.,0.) 

G=CNPLX(0.,0.) 

CC  6  1=1, M 
G=G4p*C (I) 

6  F  =  F’»ZINV 

AMP  (J)=CABS(G) 

X=PEAL(G) 

Y=AIMAG(G) 

IF  (ABS(X).GT.l.E-8)  GO  TO  60 
PHASE(J)=SIGN(90.,Y) 

GC  TO  8 

60  PHASE! J)=ATAN2(Y,X )*CCNVRT 
6  WPITE(6,20)  F( J) ,AMP(J) , PHASE! J ) ,E(J) 

20  FCFMAT(5X,F6.3,4X,E16.8,4X,E16.8,6X,E16. 8) 

10  CONTINUE 

AHAX=AMF( 51) 

CC  2C00  J=l,101 

IF  (AMP( J ) .GT.AMAX )  AMAX  =  AMP(J) 

2000  CONTINUE 

CC  2100  J=l,101 
AMPL! J)=AMP( J)/AMAX 
2100  AMPLT! J )=20.*AL0G10(AMPL( J) ) 

ITE  tl  )  =  1 
ITe(2)=0 

ECUIVALENCE(TITLE,RTB(5) ) 
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RE4L*8  TITLEdZ) 

RE4C(5,7)  TITLE 
7  FCFN4T(648) 

C4LL  CR4WP(101,F,AMPLTtITe,RTB) 

ITe(l)=3 

ITE (2)=3 

CALL  CRAWP(101,F,BA,TTB»RTB) 

STCP 

ENC 
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APPENDIX  F 


TAPPING  RESISTANCES  FOR  NON-RECURSIVE  FILTER 


out 


PIGDRE  P.1  Coefficient  Implementation 


The  contribution  of  a  single  tap,Vj^,  to  the  output, 
is  the  product  of  and  the  coefficient,  Aj^,  or 


I I  _  ^f 

I  Vk  I  ■  r  +  Rk 


The  expression  indicates  that  Rk  is  inversely  proportional 

to  Ak>  therefore  the  largest  coefficient  will  correspond  to 

the  minimum  resistance,  R„-„,  so 

min 

A  = 

max  r  +  R^in 


Solving  this  equation  for  R£  and  substituting  that  result 
into  the  Aj^  expression  yields: 

I  ^  I  _  \iax  ^  ^in^ 

'^k'  r  +  Rj^ 
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or, 

~  *  Rmin)  ■  r 

For  the  TAD- 12,  the  output  impedance  for  a  tap,  r,  is  about 
5  Kilo-ohms. 
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